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Foreword. The usual way to deal with the so-called shoch phenomenon in compress- 
ible fluids is the following. On the one hand there is the fact that in many cases of observ- 
able flow there exist narrow zones across which pressure, density and velocity undergo 
rapid changes. On the other hand, it is well-known that the differential equations of 
inviscid perfect fluids fail to supply solutions satisfying certain boundary conditions that 
can be realized physically. One therefore makes the assumption that these differential 
equations are valid in regions of the (x, y, 2, t) space which are separated from each other 
by discontinuity surfaces whose shape is a priori unknown. From physically plausible 
hypotheses one then derives necessary conditions for the values assumed by the physical 
variables on either side of the discontinuity surfaces. Such conditions were first given by 
Riemann, and later modified by Rankine and Hugoniot. It is finally assumed—and con- 
firmed at least in special cases—that the differential equations combined with these 
transition conditions are sufficient to determine both the discontinuity surfaces them- 
selves and the continuous flows in the regions between them (see [4], pp. 116-118, 134- 
138). 

A different approach, as suggested by R. von Mises to the author, is followed in the 
present paper. The sole basis is formed by the system of partial differential equations 
(Navier-Stokes equations) which govern the motion of a viscous, heat-conducting, 
compressible fluid. No additional assumptions of any kind are introduced. It is proved 
that the integrals of these equations include a class of solutions of the boundary layer 
type; that is to say, solutions which asymptotically converge (with vanishing viscosity) 
towards flow patterns entirely different from those which are obtained when from the 
start viscosity is neglected. For a small viscosity coefficient u, these flows have rapid 
changes of the physical variables across certain narrow regions, the widths of which con- 
verge to zero as u — 0. In the limit the values of the variables on the two sides of the 
transition are subject to equations which are identical with the Rankine-Hugoniot con- 
ditions. These conditions, obtained here without any hypotheses, are thus proved to be 
not only necessary but also sufficient for the existence of shocks. 

1. Introduction. In a recent paper [6], R. von Mises has discussed the occurrence of a 
shock in the one-dimensional steady flow of a reai fluid. Here “real’’ will be used to 
denote ‘“‘heat conducting, viscous and compressible.’’ If we write his fundamental equa- 
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tions (1) through (5) in a dimension-free form by referring p, u, 2, p to standard values 


3 
Po» Uo, A, pollo , We Obtain: 





pu = m, (1.1) 
mu+p—o = Cm, (1.2) 
(f 4. T \4up- 0) -H=C 1.3 
m\5 gre ulp—a = Cym, (1.3) 
p = pT, (1.4) 
with 
_ 4% 
°  8R dx’ 

(1.5) 

a ae _1 di 


where 7 is referred to u5/R (R is the gR of [6]), and we use the same symbols to denote 
the non-dimensional quantities, but change the y of Reference 6 to « for convenience. 
Here RF is the Reynolds number, 


R = Pood (1.6) 
be 
and P the Prandtl number, 
P=— er Re. ; (1.7) 


These equations are based on the assumption that the fluid is a perfect gas (1.4) with 
constant specific heats. Further, von Mises considers that P, u, k may vary. For conven- 
ience we shall also assume y, k and therefore P are constant, though the results are equally 
true without this restriction. Finally, if we write 


(1.8) 


the equations (1.5) become 





4 du « id ; 
o 3 ds’ ~ «—-1P ds (1.5') 


and the set of equations (1.1) through (1.4) and (1.5’) is now independent of R. 

In order to discuss the boundary layer nature of shock transition, we recall the as- 
sumptions of the Prandtl boundary layer theory, according to which the Navier-Stokes 
equations are reduced to a new system of “‘boundary layer’’ equations (see [2]). If we 
denote by (v; , v2) the components of the velocity along the normals and parallels to a 
boundary line S in two-dimensional steady flow, and by 0/ da, 0 /d8 differentiations along 


these directions, then symbolically 


t) 


Ba 


O(R™); aa = Ol), (1.9) 


4 = OR”); v, = O(1), (1.10) 
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and other physical quantities are O(1). Thus we stretch the normal distance by a factor 
R'’’, magnify v, , the normal component of velocity, by R'”, and then consider that all 
new quantities and new derivatives have the same order of magnitude O(1). Neglecting 
all but the terms of highest order in our equations, we obtain the reduced set of ‘boundary 
layer’”’ equations. A solution of the latter equations constitutes an asymptotic integration 
of the Navier-Stokes equations. Now consider shock transition at S. Here, in contrast to 


(1.9), (1.10), we assume 


0 0 

— = O(R); — = Ol 1.9’ 

y= O(1); v2 = O(1), (1.10’) 
where these are suggested in the following manner. If we superpose a constant velocity 
v. = O(1) on the motion discussed by von Mises, then we obtain a two-dimensional 


steady motion with S a straight line. For this motion, (1.9’) and (1.10’) clearly hold, in 
view of (1.5). Thus, in the general case, we stretch the normal distance by a factor R and 
then consider all new quantities and their new derivatives to have the same order of 
magnitude O(1). Neglecting all but the terms of highest order in our equations, we obtain 
the reduced set of ‘shock transition’? equations, a solution of which constitutes an 
asymptotic integration of the Navier-Stokes equations. These “shock transition’’ equa- 
tions are found to be completely analogous to the system (1.1) through (1.4) and (1.5’). 
Hence the existence of “shock transition’ regions in a real fluid can be inferred; further, 
the fact that as R — © these regions go over into discontinuity surfaces for which the 
Rankine-Hugoniot conditions of inviscid flow apply, can be shown. 

Having discussed the points on which the boundary layer and the shock transition are 
analogous, we must next note two important differences. First, compressibility is an 
essential for the latter, but not for the former: shocks do not occur in an incompressible 
fluid. Second, time-variation of the motion produces essential changes in shock-transition, 
but not in the boundary layer. In the former, S has its own motion; in the latter, S is 
fixed. Hence the purpose of the paper is formulated precisely as follows: Given an arbitrary 
hypersurface S (which satisfies loose regularity conditions), a V* in (¢, x, y, 2) space 
it is possible to solve asymptotically the Navier-Stokes equations so that physical vari- 
ables undergo rapid changes normal to S, and relatively slow changes parallel to S, in the 
neighborhood of S, provided that R is large enough. Moreover, in the limit R ~@, S 
becomes a discontinuity surface across which the Rankine-Hugoniot conditions are 
fulfilled. To make the discussion clearer, we shall work through two simpler cases first, 
one of steady motion and the other of unsteady. 

Finally, we shall use the fact (inherent in von Mises’ discussion) that du/ds, d’u/ds”, 
dT /ds, d’T/ds° --+ are of the same order as u, 7 --- in some fixed interval [s], depending 
only on m, C, , C2, of s, and that the end points of this interval correspond very nearly to 
the Rankine-Hugoniot conditions which strictly hold for s ==F ©. This is equivalent to 
the result that the “thickness” of the shock is extremely small in von Mises’ case, and 
that outside conditions are very nearly uniform. 

2. Steady two-dimensional flow. The momentum equation in vector form may be 


written non-dimensionally : 


]. s my 
ow x V+ grad } y ) = grad p + ap grad 3 R cule (2.1) 
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where the vorticity w and the dilatation 8 are given by 
o = curl V, ’=div V (2.2) 
({1]}true also for compressible and heat-conducting fluid). In (2.1) R is the Reynold’s 


number and the other quantities have their usual meaning. The equation of continuity 
may likewise be written: 


div pV = 0 (2.3) 

and the equation of state: 
p = pT. (2.4) 
The energy equation is not commonly derived under such general conditions but a short 


computation yields (appendix A): 


a a 
pV - grad (45 +L) 


= £1 1 om — aiy  ( pe} ov | pigy laa 
-—— iv I\?P T $5 R a 2); 


where P is the Prandtl number (assumed constant) and the dilatation # is given in (2.2) 
together with the vorticity o. 

These are the equations of motion and we wish to refer them to a particular set of 
orthogonal curvilinear coordinates having coordinate lines normal and parallel to the 


4 const 


eS 


cont 














peconsf. peconst. 


Fia. 1. 


given line S, (Fig. 1). This system of coordinates was used by von Mises [2] to obtain the 
boundary layer equations for an arbitrary curved boundary line S. The square of the line 
element in these coordinates may be written 


ds? = da? + (r+) dé”, (2.6) 


where £ is the arc length on S measured from an arbitrary fixed point, r = r(@) is the 
radius of curvature of S, and a is the normal distance from S. Hence if we write the 
components of velocity referred to these axes as (v, , v2) we have, instead of (2.1) and 
(2.2), 
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re a}. , 49 1 _4 de 
fo B+ to, 298 r+ T BR aa Rr+a0p’ (2.7) 
Ov, r Os ts a8 , 1 dw 
aaa: nee ia 
0 (r+ 
~ F : a abs + Zo}, 
(2.9) 
7 2 (rt2,,) - 2, 2 wh, 
oo? 2 a 2( r "1 
and instead of the continuity equation (2.3), 
tGts J 4 tege 
3 ( A) 9B (pv2) = 0 (2.10) 
The values of the various terms in the energy equation (2.5) may be obtained from 
+) (bterkeddr ee) 
Vv: grad (3 7 = aa tr 4+ aap , Sr (2.11a) 
2 r éo rt+e2) 2( r 2) 
¥ -2 ( r + 0B \r + a dB (3.116) 
V X o = (bw, —vw) (2.11c) 


with use of (2.9), i.e. we remember, for instance, that the expression for 3 gives us the 
divergence of a vector V. 

In agreement with the ideas of Sec. 1, we now replace a by s/R and consider in our 
new system of coordinates that all functions and their derivatives are of the same order of 
magnitude. Then pint only the highest order terms in 1/R, we may write the 
previous equations (2.7) through (2.10) as: 


a, _ _ op ,4a (2) ’ 
oe 8 + 3 ds (2.7') 
% _ 9/(w r 
m1 as a8 Oe} 
Ov; v2 


= (2.9’) 


0 ’ 
3, (err) = O (2.10’) 


and the energy equation, from (2.1la) through (2.11e), as 
a (1 ys Ej} ata. Lar _ al ( 22),, | 
me (dvre = Pas asL\? * sR” 


+) - 53) 





(2.11) 
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It is clear immediately that (2.8’) with (2.9’) is satisfied by 


v, = v2(8) (2.8”) 


which is the mathematical expression of the fact that the component of velocity parallel 
to the line S remains unaltered along a direction prependicular to S. The equations (2.7’), 
(2.10’) and (2.11’), with use of (2.9’) may be immediately integrated to give 


pv, = m8), (2.10’’) 

4 0 1 J 4? 

m, + p— 3 = = ¢,(8), (2.7°") 

=o = ( Saw) meres 
m3 tile, soe | + \P — 3 on) . Pas ~ 24): (2.11") 


where we have made use of (2.8’’). The equations (2.4), (2.7’"), (2.10’’) and (2.11’’) are 
equivalent to (1.1) through (1.4) and (1.5’), with v, for u and ¢, , c. for mC, , mC, . 
These equations have been obtained under the sole assumptions that (i) F is large, 


ty 











Fia. 2. 


and (ii) changes along the shock are of the same order as du/ds, 0°u/ds° , 0°T'/ds° , ete. 
Thus they are accurate for the interval [s] of section 1. Hence the existence of a shock is 
demonstrated. Letting R —@, the equations become exact without assumption (ii) and 
[s] may (in this case alone) be extended to 


Thus S becomes a discontinuity and the Rankine-Hugoniot conditions hold. 
3. Unsteady one-dimensional flow. The equations governing the motion are here, in 
non-dimensional form, 


re] re) . 
ay (Pu) + a (ou + p — 0) = 0, (3.1) 


#0 
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dp, Alou) _ 
at + ox a 0, (3.2) 
Sifts, 2. fis. f= ) se 
ot {ol} + on :)} + Ox {on( + toe | +Ht+(e pu} sal (3.3) 
with 
4 du «-11 oT 
a HH = 4 PR ax (3.3a) 


and the equation of state is 


p = pT. (3.4) 


Our given surface S reduces to a line in the (z, ¢) plane and we set up the same system of 
coordinates as in the previous section (Fig. 2). Thus a = const. represents a motion fixed 
relative to the “shock”’ S, and 8 = const. represents motion directly across S. Now 
equations (3.1) through (3.3a) may be written as 


a a (24 : dA). ) - if 
a SID Y 9a oO ‘+> a8 cos ¢ + a8 sing} = 0, += 1,2,3 (3.5) 
Ai = pu +p-—a, 2 = pu, 
— oa 
a; =— pu, Ay PP; (3.5a) 
3 3 1 2 T 
A; = pu(d a4 = + H+ (¢ — p)u, A, = (3 uo + —Z_), 
2 c= i 
and 





4 {me ing + ou cos gh 
ee” Fae Be? 


_ (1) 1 JaT . ee 
H = {2 sing + —— ap cos oh, 


(3.6) 





where we have used the rules for differentiation under change of variables. We note that 
¢ = ¢(8) (3.7) 


We now replace a by s/R and consider that in the new system of coordinates all func- 
tions and their derivatives are of the same order of magnitude. Then considering only the 
highest order terms in 1/R, we may write (3.5) and (3.6) in the form: 

te) ‘ ‘ . ’ 
ae (A; — A; cot g) = 0, ¢ = 1, 2, 3. (3.5’) 


and 





4a me (e-D1 oP 
o¢= g8ney., H= : PR®™ ? 35 (3.6’) 
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If the velocity cot ¢ of the shock is written as c, the equations (3.5’) reduce to 


mu+p—o=c,(p), p(u — c) = m(8), 
; r (3.5’’) 
md u? + | — H + (p —- ou = C2(B). 

The equations (3.5’’) and (3.6’) are equivalent to those of section 1 since (3.7) holds. 
Thus the discussion at the end of section 2 may be repeated. 

4. Unsteady three-dimensional flow. In this, the most general case, the situation is 
a little more complicated. The differences between the present analysis and the preceding 
will be noted at each stage. First, the momentum equations in (t, x, y, 2) space have no 
simple “‘vector” form which allows of simple transformation to new axes, as with (2.1) in 
(x, y) space. In fact, they may probably be dealt with best in the form: 


(2 + V - grad v) = —grad p + an grad 3 + 5 VV (4.1) 
with 
S=divV, V=(u»,w). (4.18) 
The equation of continuity 
Op : 
at + div pV = 0 (4.2) 
likewise has the same drawback. The equation of state 
p = pT, (4.3) 


and the energy equation (Appendix A) 


OT + aiv ay) = £11 or - aiv | (¢ 2 a)v | 
ry + div (IV) = PR VT div | p+ aR v}V 


kK 
(4.4) 
a? roe 
+ oR VV + R div (V - grad V) 
with 
Dias. sh z 
o = curl V, [= (2 V° + ft) (4.5) 


complete the system of equations governing the motion. The energy equation also is not 
of vector form. 
Suppose that the Euclidean metric 


ds’ = dt’? + dx + dy + dz (4.6) 


space. Then it is possible to set up a system of orthogonal 


is imposed on the (#, x, y, 2 
curvilinear coordinates (a, 8, y, 6) such that (4.6) reads: 


ds? = da® + (h, + Aya)” dB? + (hz + Ava)” dy? + (hs + Aga)* dd” (4.7) 


. A a on y ; r3 
where the lines 8, y, 6 = const. are normal to the given manifold S (here a V°) and the 
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3-spaces a = const. are parallel to S (see Appendix A). Moreover h, , hz , hs , Aa, Az y As y 
are functions of 6, y, 6 alone, and 


ds’ = hi de’ + hi dy +h, d& (4.8) 


is the metric of S. To illustrate the nature of this system, we note that in the two- 
dimensional case of unsteady flow the lines (8, y const.) are normal to a surface S whose 
lines of curvature are a = 0, 8 = const. and a = 0, y = const. Also in this case h,/), , 
h./d» are the principal radii of curvature of S, and are functions of (8, 7). 

We now wish to introduce a notation for vectors A ; in our 4-space. We shall denote by 
(A, , Ao , Az , As)—round brackets—the components of A; referred to the (¢, x, y, 2) 
axes, and by {A{ , Ai , Af , A4}—curled brackets—its components referred to the 


(a, B, y, 6) axes. Moreover we shall denote by div* the divergence i.e. 


OA, , Ay, OAs , OAs 
ot ¥ Ox , oy + Oz 





(4.9) 


div* A; = 


and similarly the gradient by grad*. Finally we shall suppose that the relative orientation 
of the two sets of axes is given by the scheme 





t x y Zz 
al k, l, Ms Ns 
B| ky l, M3 Ns (4.10) 
y| ky l, Mz Ns 
5| ky l, Mz M4 





in which all elements are independent of a. 
Consider now the first of the momentum equations (4.1). We have 








= + V- grad u = (I, u,v, w) - grad* u 
(4.11) 
ees ee oe oe 
oye Oa’ hy + Aya OB’ he + Aga Oy’ hz + Aza 85)’ 
where 
v», = ky + Lut mo + ny, vo = kp + lu + mv + nw, 
(4.11a) 


v, = kz + lu + mw + n3w, Ve = ky + lu + mw + nw. 


As in the previous sections we now replace a by s/R and (considering all functions and 
derivatives to be of the same order of magnitude in the new coordinates) retain only the 
lowest powers of 1/R. Thus: 

Ou Ou 


rT] + V- gradu = Ro, as (4.12) 
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Also, to the same order, 


— RI, 2 WW _ py 2%. 
an Rl, Os ’ a ith, Os (4.18) 


Further we have 
div* (0, = 0, 0) = div* {1,U, hU, LU, 4U} 


with 


Ou l, du Le Ou # ou 


jz=fI — — = 
C Oa h, + Aya 0B he + Axa dy hs + Aga 06, 


so that 


ou 1 _ ; . 
ax? ~ (hb Aw) he + Aeahtg Fb Aca) ‘2 [((hi + Ara)(he + Ava)(hg + Aga), U] 


+ s [(he + doa) (hs + A3;a)l,U] + ce + ++} 


(4.14) 
en 
=G — ' 
Os 
to the highest order. We thus have from (4.14) 
2 e, » OU - 
Vie =i (1 — &) 3 (4.15) 


Os 
To complete the discussion of (4.1) we have from (4.1a): 
v = div* (0, u, v, w) = div* {v, — hk, , v2 — ke, vs — hg Og — By} 


(4.16) 


The z momentum equation, in (a, 8, y, 6) coordinates, is thus, 


Ou , Op ] , >» OU m 
5, — = Ss + —1,—3(@, — &) + (1 — &) a3 4.17a 
pv; as 1 as oe "! ag 1 1) 1 as” ( ) 

The y and z momentum equations are obtained in a similar manner: we merely change u 

to v and w, and /, to m, and n, , respectively, to form two new equations (4.17b) and 

(4.17c). Unlike the case before, these equations are not component equations referred 

to the new system (a, 6, y, 6) but are referred to the old system in new system terms. 


The continuity equation (4.2) can be written as 
0 = div* (p, pu, pv, pw) = div* {pv, , pus, pus , pra} 


and this becomes 


< (pn) = (4.18) 
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From the preceding discussion we may list immediately the following reductions of the 


terms of the energy equation (4.4) 


ol ‘ 0 
ry] + div (IV) = R.. (Iv), (4.19a) 
wT = Ra — 2, (4.19b) 
Os 
V*(V’) = R(1 — &) - aa (4.19¢) 
. @ 
div (p + 20)y = R 2 (p + 2), —k »|. (4.19d) 


Moreover similarly to (4.12) 
ov 
V - grad V = RY, — k,) “4 


so that similarly to (4.16) 


° »2 Cy) a a — 
div (V - grad V) = Ri a5 (6. k,) 3s (vy, ky} 


(4.19e) 
— p92 ily 25 
=R as ‘3 (v, — k,) \ 
Using (4.19a) through (4.19e) and (4.16), we may replace (4.4) by 
af (2 ‘ r\)-s21! 7 =? i tal 
Os [™ 2 eo = 3 7 tc (1 hi) Os [pe ky) 
(4.20) 
1 » A(V*) 1 a 
+ha— m+ 126, — Ky. 


Before discussing these equations, we must introduce simplifications. If we write 
2 2 2 2 2 
RPr=h+m+272,=1-& 


t is clear that n = (1,/D, m,/D, n,/D) is the unit normal in (2, y, z) space to the instan- 
taneous position of the ‘‘shock.”” Moreover —k,/D is the normal velocity of the “shock” 
in this space. Denoting then the normal component of velocity by u, and the normal 
velocity of the “shock” by U i.e. 


lLu+ mv + nw = Du,, —k, = DU. (4.21) 


We may rewrite (4.17a), (4.17b), (4.17c¢), (4.18) and (4.20) in the more familiar form: 


ou _ -h pil, eu d*u , 
pm — Ua =D ag +34 ae +? ay? (4.17’a) 
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™ — yy = amo, 1, Tu dv a 
plu, U) =. so + 3 m4 a2 + D ae? ? (4.17’b) 
4 yo mm Ol aut, ate o 
o(u ) ea — 3% 52 + D ae? (4.17’c) 
“ ° 
5 [e(un — U)] = 0 (4.18’) 
os 
Q n(2 , T )| k-1D8T Apu) 
oo yeu — sh 2 oe, os Se 
Os Ee C N\5 I " gy x ¥F oF Os 


(4.20’) 
Da(V*) , Dd&,., 
+ 2 os” + 6 ds" (u;) 


If we multiply (4.17’a) by 1, , (4.17’b) by m, , and (4.17’c) by n, and add we obtain 


ll ) Iuln ow op 4 2 at Un) (4.21) 
Os Os >] Os 
which may be integrated by virtue of (4.18’), i.e. 
p(u, — U) = m(8, x, 4), (4.22) 
to give 
mu, + p—o = ¢,(B, y, 6) (4.23) 
with 
ad 
¢= : I a (4.23a) 


Furthermore if we multiply (4.17’a) by m, , (4.17’b) by J, and subtract we have 
) a” 
m— (mu — lv) = D=s (mu — lw 
as v) as? 1 1) 
which is satisfied if 
1 , , 
D (mu — lv) = f(6, y, 4). (4.24) 


Similarly the other components of V X n are independent of a. Hence the component of V 
tangential to the instantaneous position of the shock in the (2, y, 2) is independent of a. 
Thus we may write 

te) 


ag (V) = 


: 4.25 
as (Un) (4.25) 


2 |e 


and the energy equation (4.20’) can be written in integrated form: 


m(3 Un + —t | — H+ (p — o)u, = ¢,(8, 7, 4) (4.26) 
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with 


k—-1D0T 
H= aa P Os (4.26a) 
The equations (4.3), (4.22), (4.23), (4.23a), (4.26) and (4.26a) are equivalent to those of 
section 1 with wu, for u, and definition of m(@, y, 5) slightly changed, since D = D(, y, 4) 
from (4.10). Hence discussion at end of section 2 applies equally well here. 

Summary. Given an arbitrary hypersurface S in (¢, x, y, z) space and an arbitrary 
distribution (within the loose restriction of the constant c of von Mises’ paper) of V, T’, p 
and derivatives on it, we have the following: 

(¢) It is possible to construct a family of solutions of the Navier-Stokes equations 
having large normal derivatives on S and (comparatively) small derivatives parallel to S, 
provided u is small enough. The solution is analogous to a solution of the steady one- 
dimensional equations of motion, in which the z-differentiation is replaced by 


9 _ 72 

an At’ 
where U is the normal velocity of the ‘“‘instantaneous”’ section of S by (a, y, 2) space and 
d/dn denotes differentiation along the normal to this section. Moreover the Rankine- 
Hugoniot conditions are approximated very nearly, and the component of V in the 
tangent plane to the “instantaneous” section of S by (a, y, 2) space is conserved along the 
normal direction of S. 

As u» — O the surface S becomes a discontinuity surface and the Rankine- 

Hugoniot conditions hold exactly. 

The purpose of this paper, as stated at the beginning, was to prove the existence of 
shocks, using as sole basis the Navier-Stokes equations for the motion of a perfect fluid. 
As a by-product of the analysis, we are given a method of investigating the structure of a 
shock in the most general motion of a perfect fluid. However, if such an application is 
made, it must be borne in mind that it will only be valid when the Navier-Stokes equa- 
tions hold, that is when the fluid can be considered as a continuum, in other words, for 


wea k she cks. 
APPENDIX A 


If we denote an element of the stress-tensor of a viscous fluid by o;; and an element of 
the rate of strain tensor by e;; , the generalized concept of viscosity leads to [3] 


C3; = —(» + 2, fs ud) 6; + Meg; (A.1) 
with 
Ou Ou 1 . 
7. = —— a o = —e,, = div V. A.2 
- OX; Ox; ’ 2 (A.2) 
Here we use dimensional quantities, with yw the coefficient of viscosity and 
(x, y, 2) C1 Lo, Xq), (U, v, W) = (Uy, U2, Us). Since the total energy [4] is 


Ja (3 v2 + Ar), (A.3) 


c—-l 
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the energy equation is easily seen to be 


D (? “a RT ) a O(u;0;;) 
” a = K 28 | oe fi * 
Pri No Vv" 4 a kV°T + arr, (A.4) 
From (A.1) and (A.2) it follows, however, that 
2 \ 0 ] *) ou 
uo — F “} 4 bs 5 
(p + 3 Me ju + ar, (2 } + UU, a. (A.5) 


Inserting this into the right hand side of (A.4) and using the equation of continuity (4.2) 
on the left hand side, we obtain the dimensional form of (4.4). 


For one-dimensional unsteady flow, we have V = (u, 0, 0) and 


2 ; — —— 9 
— 3H div (8V) + 3 VV’ + udiv (V - grad V) = : m = (. au), (A.6) 


and (4.2) reduces to (3.3). 


For steady two dimensional flow 0/d/ = 0 and the transformations of vector analysis 


yield 
— uV’°V*? + uwdiv (V - grad V) = udiv (w X V). (A.7) 


Hence (A.7) reduces (4.2) to the form (2.5). 


APPENDIX B 
Let r = r(v; , v> , v3) be the position vector of the manifold S (V*) in the (¢, 2, y, 2)- 
space, where (»; , v2 , v3) form <n orthogonal system of coordinates, 1.e. 
Fee8; = [Te | | Ps | Oss (B.1) 
with 


or oi : 
—_—_ > (1 = is 2, 3) 
Ov 


i 
i 


The unit normal n to S is uniquely determined by the equations: 

n-r,; = 0, (B.2) 

n-n =1. (B.3) 
From (B.3) we have by differentiation: 

n-n; = 0, (B.4) 
and hence we may write 


r; ~ 
ra I; (B.5) 
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with suitable n;; = 7;;(v1 , v2 , v3). Moreover from (B.2) we have 
n;-r; + n-r,; = 0, 
whence 
n;-r; = n;°F; (B.6) 
or substituting (B.5) and noting (B.1) 
Nis = Ni» (B.7) 


Now suppose we have a unit vector in S defined by 


s= > 3-5 (B.8) 


t=1 | ri | 


Da, 7 = > Dame. i (B.9) 
i=1 | oa i=] j=1 i 

and this lies along s if 
>> 0:8; = ks; (j = 1, 2, 3) (B.10) 


i=l 
for some k. From the theory of matrices it is known [5] that since n,; is symmetric (B.7) 
there exists at least one set of mutually orthogonal real unit eigenvectors s“’(A = 1, 2, 3) 
° » . . . rh ° . . 
satisfying (B.10) with corresponding real eigenvalues k = k*’. Taking these directions 
s’” as the coordinate directions of a new set of orthogonal coordinates (8, y, 6) in S we 
have 


Ry = pike ; Bn, = pats » is * pli (B.11) 


with pw, , 2 , Ms functions of (8, y, 4). 
Now let n complete a system of four orthogonal directions in 4-space and @ be distance 
measured along it. Then if R is the position vector in this space 


R=r+an 
with r, n independent of a. Hence, using (B.11), we have 
R, =a, Rs = fs + ang = fz (1 + wa), 
(B.12) 
R, = r, (1 + ua), R; = r; (1 + wa). 


Hence (a, 8, y, 6) form a system of orthogonal curvilinear coordinates with 


ds? = do? + (hy + dra)” dB? + (hz + Axa)’ dy? + (hs + \;a)° dé’, (B.13) 


where h? = r3 ete., independent of a, and \,; = h,u, ete., independent of a. 
We have thus demonstrated the validity of (4.7) and (4.8) as well as the statements 


about these equations. We have constructed V*-spaces a = const. parallel to S (in the 
obvious sense of the term) and used them to form an orthogonal system of coordinates. 
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ON DOUBLE-PULSE STABILITY CRITERIA WITH DAMPING* 


BY 
M. MORDUCHOW anp L. GALOWIN 
Polytechnic Institute of Brooklyn 


Summary. A simple necessary condition for the stability of a linear dynamic system 
with elastie and damping characteristics which vary periodically with the same period 
is derived. General explicit necessary and sufficient conditions for stability are then 


developed for a double-pulse system. Such a system can be characterized by a pair of 
eigenvalues, or complex frequencies, corresponding to each half-period, and the stability 


of this system depends only on these complex frequencies. It is shown that a necessary, 
though not sufficient, condition for the stability of any such system is that the arith- 
metic mean of the real parts of all four of the complex frequencies over an entire period 
be negative or zero. This is shown to be true, more generally, for an n-pulse system. 
The physical significance of the results is discussed, and numerical examples are given. 

Introduction. \ variety of physical phenomena’ are characterized by a linear homo- 
geneous differential equation with periodic coefficients. Sometimes the chief interest is 
in the stability of the solutions of such equations, rather than in the solutions themselves. 
If the system which these equations characterize is such that the coefficients vary 
periodically) by changing discontinuously from one constant value to another one or 


F(t) 








| 
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Fic. 1. Double-Pulse Function F(t) with Period T (= 27). 


more times during a period of variation (Figs. 1 and 2) then an exact solution of the 
equations can in principle be obtained without difficulty. Coefficients varying in the 
manner of Figs. 1 and 2 are called “rectangular pulses’’. In particular, if the coefficient 

*Received March 29, 1951. 

1Examples are given in: (a) B. Van der Pol & M. J. O. Strutt, On the stability of the solutions of 
Mathieu’s equation, Phil. Mag., 5, 18 (1928); (b) S. Timoshenko, Vibration problems in engineering, 
second edition, D. Van Nostrand Co., N. Y., 1937, Chapter III; (c) J. P. Den Hartog, Mechanical vibra- 
tions, third edition, McGraw-Hill Book Co., N. Y., 1947, p. 408. 
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assumes only two different constant values, one for each half-period of variation (Fig. 1) 
it is called a ‘‘double-pulse”’. In cases in which the periodic coefficients do not vary as 
rectangular pulses, it is often possible to obtain approximate or at least qualitative 
solutions by replacing these coefficients by either double pulses or four pulses.” This 
type of approximation has, in fact, been recently applied.” 

The purpose of this paper is to present explicit stability criteria for systems which 
can be represented by equations with periodic, especially double-pulse, coefficients for 


F(t) 
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Fic. 2. Four-Pulse Function F(t) with Period 7’. 


the damping as well as for the elastic characteristics. Such a system without damping 
has been analyzed in detail by Meissner* and by Van der Pol and Strutt.” Schwerin” 
subsequently extended this investigation to systems with linearly varying character- 
istics. A particular double-pulse system with damping has also been investigated,’ 
without, however, attempts to derive explicit equations and conclusions for damped 
systems in general. 

The results obtained in the present analysis serve to present useful theorems of 
physical interest, to demonstrate some significant conclusions which might not be 
readily foreseen without the mathematical analysis, and to facilitate any actual stability 
calculations based on the double-pulse method. 

General theory. The free motion of a linear dynamic system with spring stiffness k 
and damping coefficient f varying periodically with time, with period 27, can be repre- 
sented by the equation: 


a’ + fie’ + k(x = 0 (1) 


2B. Van der Pol and M. J. O. Strutt, op. cit. 

3G. Horvay and S. W. Yuan, Stability of rotor blade flapping motion when the hinges are tilted. General- 
ization of the “rectangular ripple’ method of solution, J. Aero. Sci., 14, 583-593, (1947). Also, S. W. Yuan 
and M. Morduchow, On the stability of the transient motion of helicopter blades in flapping and lagging, 
Reissner Anniv. Vol., J. W. Edwards, 1949, p. 163. 

‘E. Meissner, Uber Schiittelschwingungen in Systemen periodisch verdnderlicher Elastizitdt, Schweizer. 
Bauzeitung, 1918, p. 95 (1918). 

Sop. cit. 


8. Schwerin, Ein allgemeines Integrationsverfahren fiir quasiharmonische Schwingungsvorgdnge, 


Z. Techn. Phys. 12, 104 (1931). 
7G. Horvay and 8. W. Yuan, op. cit. Also, S. W. Yuan and M. Morduchow, op. cit. 
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where x(/) denotes displacement, the primes denote differentiation with respect to the 
time ¢, and f(t) and k(¢) are functions having the period 27. By letting 


x = v(t) exp (- / fat), (2) 


Eq. (1) can be transformed into the form 

v’ + k,(tv = 0, (3) 
where k,(t) has the period 2x. By Floquet’s theory” at least one particular integral of 
Eq. (3) is of the form 
v, = o(t) exp (ut) (4a) 
where ¢(/) is a function with the period 27, and u is a definite constant. From this it 
can be shown that a second independent integral of Eq. (3) must have the form 

v. = ¥(t) exp (—ul) (4b) 
where y¥(/) has the period 27. 

From Eq. (2), and from the fact that, due to the periodicity of f(d, 


pt+2r 


| f(t) dt = constant = C, (5) 


say, it follows that the solutions x,(0) and 2,(t) of Eq. (1), corresponding to v, and 0, 
respectively, will satisfy the following relations, regardless of the time ¢: 


x(t + 2r)/z,(t) = 11, x(t + 2Qr)/2x.(t) = oe (6a) 
where o, and g, are constants given by 
o, = exp (2ru — 30), o. = exp (—2au — 3C). (6b) 
The general solution of Eq. (1) can be written as 
x = A,z,(t) + A22,(2), (7) 


where A, and A, are arbitrary constants. 

The dynamic system may be defined as stable when both of its component free 
motions x, and x, subside with time, and unstable if either of these motions amplifies 
with time. From Eqs. (6a) it then follows that the system will be stable if | o,| < 1 
and | o,| < 1, while the system will be unstable if either | o, | > lor|o.,| > 1. The 
o,| = land | o,.| = 1 may be defined as one of neutral stability. 

Equations (6b) imply that o,¢. = exp (—C). Hence in order that both | ¢,| < 1 
< 1, it is necessary that C > 0. This can be expressed by the following theorem: 


case 


and | op 
TuHEeoreM I: A necessary condition for the stability of a dynamic system characterized by 


an equation of the form (1) is that 


[ soazro 


8M. G. Floquet, Sur les équations différentielles linéaires a coefficients périodiques, Ann. Ec. Norm. 12, 
47 (1883). Also, E. T. Whittaker and G. N. Watson, Modern analysis, Macmillan Co., N. Y., 1948, pp. 


412-413. 
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Physically this condition can be interpreted as requiring that the average value of 
the periodic damping coefficient f(t) during any period must be non-negative. Thus, for 
stability, there must in the average be a positive damping force. 

Double-pulse system. To investigate, in particular, the stability of a double-pulse 
system with damping, let f(t) and k(t) vary periodically in the manner of Fig. 1, i.e. 
let f = f, andk = k, for 0 < t <z, while f = f, and k = k, for r < ¢ < 2m; the 
constants f, , f2 , k, and k, may be positive or negative. Then Eq. (1) can be solved 
for each of the two time-intervals as a linear equation with constant coefficients. The 


solution can be expressed as: 


x = C, exp (p,0) + fi exp (pt), OStsr 
(8) 


x = C, exp (pot) + Cz exp (put), xs t<2r 


where the C; and C; are arbitrary constants, the p; and p; are constants given by: 


The p; and 9; are essentially eigenvalues corresponding to each half-period of variation 
of the damping and elastic coefficients, and are sometimes called ‘‘complex frequencies”. 
Their physical interpretation follows from Eqs. (8). Thus, if a system were characterized 
by constant dynamical coefficients leading to a (single) pair of complex conjugate values 
for p and @, then the system would oscillate freely with a natural frequency proportional 
to the imaginary part of p or p, and would be damped out with a logarithmic decrement 
directly proportional to the negative of the real part and inversely proportional to the 
imaginary part of p or p. A system with constant coefficients characterized by a single 
pair of real values of p and p would be non-oscillatory and would be stable if both p 
and 7 are negative, but unstable if either p or D is positive. For convenience, the negatives 
of the real parts of p; and 9, will here be called ‘damping factors’”’. 

The arbitrary constants C; and C; in Eqs. (8) are related by the two conditions of 
continuity of displacement x and velocity x at t = 7. Moreover, the following two 
additional conditions can be prescribed: 


(2), mor = o(Z):00; (2%’) <3 = o(2’),<0 ; (10) 


where o is a constant, to be determined. It can be shown from Eqs. (2), (4a, b) and (7) 
that by satisfying conditions (10), initial conditions are implicitly prescribed so that 
A, or A, in Eq. (7) will be zero,’ and so that consequently the two values of ¢ which 
will thus be obtained will correspond exactly to o, and o, respectively of Eqs. (6b). 
By applying the continuity conditions and conditions (10) to Eqs. (8) a set of four 
linear homogeneous equations in C, , C, , C. and C, is obtained. The condition for the 
existence of a non-trivial solution to these equations leads to a determinantal equation, 
which when expanded can be written in the form: 


o +No+M =0, (11) 
where 


*This fact, although implied, often appears not to be clearly stated in the literature. 
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_ Ps = B)(p2 — p)(PsPs + PrPs) + Po — PIG» — P)P:P2 + PsP2) 99) 
(Pi — Di)(D2 — Po) 





N 


M = P,P,P.P, , (12b) 


P; = exp (p;m), P; = exp (p;n). 


r r\2 1/2 
c= -> + | (¥) - au | (13) 


of Eq. (11) then determine the stability of the system. 

Three general cases of double-pulse systems, depending on whether the eigenvalues 
p; and p; are regl or complex, can be distinguished. By substitution into Eqs. (12a, b) 
the following expressions for N and M in these various cases are obtained (a; , }; , ¢; 
and d; denote real positive or negative constants, while i = (—1)’”). 


The roots 


Cass I. Two pairs of complex conjugate eigenvalues: 


Dp, = —a, + Dt, Di = —a, — 51, Po = —a, + byt, D2 = —a, — byt. 
2b,b.N = (exp [—72(a, + ae)]){[(a, — a2)” + (b; — b2)] cos (b, — b.)r 
(14a) 
= ii = 2)” + (b) + b.)”] cos (b; + b2)r} 
M = exp [—2z(a, + a,)] (14b) 
Case II. Two pairs of real eigenvalues: 
“= wy Pi -_ —d, ] P2 = —C2 , Do _ —d, , 
(c, — di)(c2 — d.)N = (d, — d.)(c. — ¢,) {exp [—7@, + d,)] 
+ exp [—z(c, + d,)]} (15a) 
+ (co — d,)(d2 — ¢,) {exp [—r(, + c2)] + exp [—x(d, + d,)]} 
M = exp [—7r@ + d, +c + d,)] (15b) 


Case III. One pair of complex conjugate, and one pair of real eigenvalues: 


Di = —a, + by, Pi = —a, — byt, P2 = 2, pP2 = —d, 
b,(d, — c.)N 
= (exp [—7a,]){[(d2 — a:)(a1 — ¢2) — bi][exp (—xd2) — exp (—me,)]sin br — (16a) 
+ b,(c. — d2)[exp (—md,) + exp (—2c2)] cos br} 


M = exp [—7(2a, + c. + d,)] (16b) 
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For any given double-pulse system, the stability characteristics can be quickly and 
straight-forwardly determined by calculating the p; and p; from Eqs. (9), evaluating 
N and M by means of the appropriate pair of expressions from the set (14a, b), (15a, b) 
and (16a, b), and then obtaining the two values of o from Eq. (13). It should be noted 
that these calculations will always involve only real quantities. 

From Theorem I it follows that for stability of a double-pulse system it is necessary 
that (f, + f.) > 0. According to Eqs. (9), this is equivalent to the following theorem: 
THEOREM II: Jf a double-pulse system is characterized in one half-period by eigenvalues 
p, and p, and in the other half-period by eigenvalues p, and p, , then a necessary condition 
that the system be stable is that the arithmetic mean of the real parts of p, , P, , Pp» and pz 
be negative or zero. Alternatively, a sufficient condition that the system be unstable is that 
the arithmetic mean of the real parts of p, , Pi , P2 and py. be positive. 


By virtue of Theorem I and of relations of the form (9) it can be quite similarly 
proven that Theorem II is valid more generally for an n-pulse system. Moreover, it 
similarly follows that if during any period of an n-pulse system, f(t) assumes constant 
values f; , fo, --: , f, during unequal time-intervals, then a necessary condition for 
stability is that the weighted mean of the real parts of the complex frequencies Pi to’ 
p, be negative or zero, the weights being proportional to the time-intervals corre- 
sponding to each of the p’s. 

It may be noted that the above necessary condition for the stability of a double- 
pulse system is in a sense less stringent than that of a system with constant elastic 
and damping coefficients, since in the latter case it would be required that each damping 
factor be positive, and not merely the average of the damping factors. It must be 
observed, on the other hand, that the theorems derived here give only necessary, but 
not sufficient, conditions for stability. This follows essentially from the fact that although 
the conditions | ¢, | < 1, | o2 | < 1 always imply o,0, < 1, the converse is not necessarily 
true. As a consequence, a double-pulse system may be unstable even under conditions 
for which a constant-coefficient system would be stable. For example, it is well known 
that in the absence of damping forces a constant-coefficient system will be neutrally 
stable when the restoring force is positive (k > 0), leading to purely imaginary complex 


frequencies. However, as can be verified by putting a; = a, = 0 in Eqs. (14a, b) and 
noting that values of b, and b, (e.g. b; = 1.5, b. = 0.5) exist for which | VN | > 2 and 


> 1, a double-pulse system may in such a case be unstable. 


hence | o¢ 
14) and (16) in conjunction with (13) it is seen that although the neces- 


From Eqs. 
sary condition for stability stated here involves only the damping factors, the necessary 
and sufficient condition involves the natural frequencies as well. This is partly in contrast 
to a system with constant damping and spring characteristics, where the stability de- 
pends exclusively on the damping factors. Considering, for example, Case I, which is 
in practice the most important case, it can be seen from Eqs. (14a, b) that a large positive 
value of (a, + a.) tends to keep the values of | N | and | M | far below unity, and hence 
strongly tends to make the system stable. This result appears quite reasonable on 
physical grounds. Nevertheless, it will be found that for any given value of (a, + a2), 
i.e. for any given amount of damping in the average, there will exist values of the fre- 
The theorem which follows could also be 


0 » directly derived by noting that for stability it is necessary 
that | | < 1, and by using the expressions 


(14b), (15b) and (16b) for M. 














ee 


ce 
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quency ratio 6,/b, for which the system will be unstable. This can be verified by ob- 
serving that for any given positive value of (a, + a,), however large, sufficiently large 
values of b,/b. (or 6./b,;) can be found for which | N | as given by Eq. (14a) will be 
large enough to make at least one of the roots in Eq. (13) greater than unity in absolute 





value. 
EXAMPLES 

1. As an example illustrating the insufficiency of the necessary condition derived 
here for stability, suppose that in Case I a, = 5, a. = —4 and b, = b, = 1.5. Then 
the average of the damping factors a, and a, is here positive; nevertheless this system 
is unstable, for Eqs. (14a) and (14b) yield: N = 38.0 exp (—7), M = exp (—2z). The 
larger of the roots of Eq. (13) will be |o| = 1.64, which exceeds 1 and hence implies 
instability. 

2. The purpose of the second example is to show that a system pertaining to Case IT 
can, contrary to systems with constant characteristics, be stable even when the average 
of a single pair of damping factors over a period, e.g. (c, + d.)/2 or (e,; + ¢2)/2, is 
— 1.333, d, = 1.333 and d, = 0.880. Then (c, + c.) < 0. 


negative. Thus, let c, = 1, co = 
Eqs. (15a, b) however give: N = —0.209, M = 0.00169. The roots (13) for ¢ are o;, 
0.201 and o 0.0081, and the system is therefore stable. 


It may be noted that in either of the above two numerical examples, it might have 
been difficult to predict in advance the stability or instability of the given system. 


CONCLUSIONS 
1. A necessary condition for the stability of a linear dynamic system with damping 
coefficient f(4) and elastic coefficient both varying periodically with the period 27 is 


2 


[- f(j) dt > 0. 


Jo 

2. In particular, a necessary condition for the stability of an n-pulse system is that 
the weighted arithmetic mean of the real parts of the 2n “complex frequencies” char- 
acterizing the system during any period of the variation of the damping and elastic 
coefficients be negative or zero. The weights must be proportional to the time-intervals, 
of any period, corresponding to each of the complex frequencies. For equal time-intervals, 
of course, the simple arithmetic mean is taken. The case n = 2 with equal time-intervals 
(‘“‘double-pulse’’), was of especial concern here. 

3. Contrary to systems with constant damping and elastic coefficients, the necessary 
and sufficient condition for the stability of a double-pulse system depends, in general, 
not only on the damping factors, or (negative) real parts of the complex frequencies, 
but also on the natural frequencies, or imaginary parts of the complex frequencies. 
Moreover, contrary to a system with constant coefficients, a double-pulse system can 
be stable even when not each of the damping factors is positive or zero. 

4. The stability characteristics of any given double-pulse system with (positive 
and/or negative) damping can in practice be quickly calculated from the equations 


developed here. 
The authors hereby express their thanks to Professor R. M. Foster for his valuable 


suggestions and comments. 




















THE RANDOM VIBRATIONS OF A STRING* 


BY 
M. R. SPIEGEL 


Rensselaer Polytechnic Institute 


1. Introduction. In the theory of linear electrical networks one encounters systems of 


equations having the form 


Ling! + Regi + Grn) = DEx, goijerryn (1) 
k=1 


where the q, represent charges, L;, inductances, Rj, resistances, G;, reciprocals of 
capacitance, the /;, are random e.m.f’s**, and primes denote differentiations with 
respect to time 

The theory of such systems of equations, quite carefully examined during the war 
years in applications to noise in electrical networks, turns out also to be applicable to 
the various mechanical systems—in particular, it is applicable to the system of a vibrating 
string with fixed end points. A general theory for the system of equations (1) has been 
developed by Uhlenbeck and Wang [2]. Some of the results which we shall obtain in 
this article have been derived without making use of the general theory. We shall make 
comparison with these results and in addition we shall derive several more results. 
The results of Uhlenbeck and certain of his co-authors [1, 5], have been derived directly 


t 


from the differential equation of motion for the string. In order to apply the general 
theory to the vibrating string, it is necessary to discretize the string. We therefore 
assume that the stringt is made up of n + 2 particles, (2 fixed, n vibrating) of equal 
mass m harmonically bound together by means of massless springs. Furthermore let us 
assume that this system of particles has random forces acting on it and that as a result 
the system vibrates, the vibration taking place in a plane. As a last assumption we 
suppose that the vibration takes place in a viscous medium so that each of the particles 
undergoes a damped vibration. 

When we have obtained our results for the discretized system, we can derive the 
results for the case of a continuous string by a limiting procedure, namely that of letting 
n the number of vibrating particles go to infinity while the total mass and length of 
the system remains constant. 

In this article we shall derive the following: 

1) the characteristic function of the sum of the squares of the deviations of the dis- 
placement of the particles from their given initial positions—also the corresponding 
characteristic function for the continuous string, 


*Received April 23, 1951. This article represents part of a thesis submitted in partial fulfillment of 
the requirements for the degree of Doctor of Philosophy at Cornell University. The author wishes to 
acknowledge here his indebtedness to Professor Mark Kac of Cornell University who supervised the work. 

**The word random used here in a heuristic sense obviously needs more careful definition and 
interpretation. A more precise investigation than can be given here is found in an article by Uhlenbeck 
and Wang (see [2]). 

+We adhere to the assumptions made in the usual derivation of the differential equation of the string, 


i.e. small vibrations, constant tension throughout, ete. 
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2) the mean square deviation of the displacements of the n particles from some given 
initial distribution and the corresponding mean square deviation for the string, 

3) the mean square deviation of a single particle of the system given that initially 
it had a certain displacement 6, and the corresponding limiting case. 


2. Certain preliminary considerations. Let the displacement of the k-th particle of 
the system be denoted by y, . Let F,(¢) denote the random force acting on the k-th 
particle. By a simple application of Newton’s law of mechanics we have the result 


47 i r ; T ¢ 7 
my.’ + Kyi — 7 (Yer — 2ye + Ye-1) = F,(0, (2) 


where K is a quantity dependent on the damping and r is the tension in the springs 
connecting the particles. The tension is assumed to be constant at all times. The quantity 
d is the distance between particles, also assumed constant. 

In order to prepare ourselves for the limiting case of the string, it will be necessary 


to define a density u in such a way that ud = m. Assuming the case of n + 2 particles 
(n vibrating, 2 fixed) it is easy to see that the following relations will hold 
L M 
C= ——— C= (3 
n+ 1 n+ 1’ ) 
_m _ M 
# d L ? 


where J is the total mass of the system and L is the total length for the limiting case 
of the string. Both M and L are given constants. Also, let K = Bd where B is a constant 
dependent on the magnitude of the damping which is dependent on the viscosity of the 
medium in which vibration is taking place. 

With the use of these relations (2) becomes 


udy;’ + Bdy; — aC +1 — 2ye + yx-1) = F,(t). (4) 
¢ 


This equation is of the form (1) if the following relations hold: 
Liz = Sj,pd, R;, = 6;,Bd, 
(5) 
GQ = 2r G, se, G; i-) = —r’, 


where 6,, is the Kronecker delta, equal to one if 7 = k and zero otherwise, and r° = 7/d. 
From physical considerations it is evident that the quantity B = 8/u where 8 is 
the actual damping coefficient for the medium in which vibration is taking place. 
As Uhlenbeck and Wang [1] have pointed out in their discussion of the system (1) 
the 2n variables [9,(é), g(t), «++ , @a(Q); a1(0, 92(0, --- , g4(O] or in our case the collection 
[yi(t), --- , yn; yO, , y,(t)] will form a 2n-dimensional Markoff process governed 


by the Fokker-Planck or generalized diffusion equation 


oP <9 ) Le _¢ 
at ain > dx; (AP) + 2 2, Ox; OX, nF), " 


where P is the probability density associated with the distribution of the quantities 
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[%, , 2, *** , Lop] Here 2, , 22, -+* , 2, are the 2n variables describing the collection 
[yi(), «+ , yal); yi(O, «++ , yi(O]. As derived in the article 


A; = > AjnLe 5 (7) 


where the a;, are the elements obtained from the matrix 
0 I 
A= : (8) 
“Eg -L2]! 
and the D,, are obtained from the matrix* 
0 0 


D= (9) 
0 2«TL RL | 





In the matrices (8) and (9), L, R and G are themselves matrices consisting of the 
elements L;, . R;, and G;, first encountered in the system of equations (1) and which 
for our purposes may be defined by (5). 

In order to solve the Fokker-Planck equation it turns out to be convenient to make 


the linear transformation 


2; = +B Cir Ly > (10) 


k=1 
and it then follows that (6) becomes 
oP — 0 leo a’P 
— = d 2,P = tn 11 
at > ” 02; GP) + 2 >») ™ Oz; 2,’ (11) 
where the o;, are obtained from the matrix CDC™, C being the matrix which diagonalizes 
the matrix A and which thereby yields the eigenvalues A; (j = 1, 2, --- , 2m). 


It is easy to show that the eigenvalues satisfy the determinantal equation 


Q P 0 ove 0 0 


r Q P ore 0 0 


SF @ 0 0 
= 0, (12) 








0 0 0 cee i Q 


By 
Be 


2r” r 
ud’ sl pd 


If we denote the determinant (12) by A, we may show that the following difference 


Q=vN+—+4 


equation is sat isfied 


*Here «x is Boltzmann’s constant and 7’ is the temperature of the medium in which the vibration is 


taking place. 
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A, _ QA,-1 — P*A,g-2 (13) 
subject to the boundary conditions Ay = 1, A, = Q. Solving the difference equation 
by the usual methods we obtain 


_ 3 Q+(Q°-—4P)"\"" _ fa — (Q’ — 4P)""\""" 
A= gre || 2 j 2 | ow 


_ 


yielding the evaluation of the determinant (12). 





Upon finding the conditions under which A, = 0, we obtain the equation 

 _a.rF mr - 
Vv+—+ = — cos ———, (15) 

be ud pd » l 

and so the 2n eigenvalues \; are given by 

B | B’ Ar i ae i. 2 

Sn a 5-36 + loa =; m=1, +--+ ,n (16) 
2u | 4u pw 2n + 2 


The evaluation of the matrix C which diagonalizes A is a little too long to present here 
but it is straight-forward. It turns out that the matrix C is given by* 


\ 
. kmr ewe oe om 
Cos? “ia + Oe .  * 
: n-<r 1 | 
( m= 1, +++, (17) 
kmr 
Cie oes ae sin —— | 
sa n+ 1 } Dimy = 1, -** 5 Qn. 


After obtaining the matrix C, a simple evaluation of CDC™, where D is given in (9), 
yields the o;, of the Fokker-Planck equation (11). This evaluation gives 


2«TB(< ] k 

2xT | . mr . mkr 

i, = | 1 sin —~ sin - : (18) 
pd jae n+ 1 n+ 1/’ 

where o;, = o7+1, if j is odd and o;, = oj, if k is odd. Evaluation of the sum in (18) 


shows that the o;, are given by the elements in 


|| ae a 0 0 tee 0 











] a a 0 0 . 0 | 
| 
0 0 a a 0 | 
| 
(19) 
| 
| 0 0 a a 
| 0 0 a a 
where a = «7TB(n + 1)?/uL. 


*The significance of the bracket in p,m) is that corres yonding to the value 1 for m we have p,,,..) equal 
P(m) I ) ] 
to 1 and also 2, corresponding to m = 2 we have p,m) equal to 3 and also 4 etc. 
? I F } j 
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The transformed Fokker-Planck equation (11) has for its solution a 2n-dimensional 
Gaussian distribution with the average values 


Z, = bjo exp (Ajf) (20) 
and covariance 
(2; — Z,)(@e — %))average = ~—-4— [1 — exp {A; + 
; — 2) Ge — %1) averse = Fay. HL — exp (0s + M4] (21) 


where 6;) are the initial values of the z; . 

3. The characteristic function. Using the results obtained thus far we may now 
write the probability density function for the system of the n vibrating particles using 
a well-known result from the theory of multivariate Gaussian distributions. The prob- 


ability density has the form 
l - - 
P(z, , Ze, *°* » Zn) = R, exp {-4 .¥ By (2; — Z;)(@ — a0}, (22) 


where the B,, and FR, are constants determined from (20) and (21), and the values of 


\ are determined from (16). 
We are now ready to obtain the characteristic function of the sum 


1 : ; 
17 ate <a, (23) 


i.e. the characteristic function of the normalized* sum of the squares of the deviations 
of the displacements from some arbitrary initial values. This characteristic function, 
once we have obtained it, will at least in principle, enable us to obtain the characteristic 
function of 

1 al ‘ 

i} (X@, 0 — S@P de (24) 
which is the corresponding limiting case of the finite sum (23). Once we have this we 
are theoretically able, by means of the Fourier Inversion Theorem of Probability, to 
obtain 


PAI [ [X(2z, t) — S(x)]’ dx < } (25) 


that is the probability that the mean square deviation of the string from a given initial 
distribution be less than some given quantity e. 
The required characteristic function is defined by 


Fie) = [ see [ R,, exp EE 7. (x; — s;)° 
| (26) 


1 = - 
ee 2 > B;,(2; ‘aie Z;)(2 _ 2) | dV 


*This normalization is necessary if our results for the limiting case of the continuous string are to be 


meaningful. 
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where the integration element dV = dz - dz, . If the integration in (26) is performed, 


we obtain the result 


where 
] = mkn 
U, = —— > &G,, — s,) sin ——, (28) 
a-—l &S n-+ 1 
and 
V. = [1 — Pauus.al” (29) 
Here, for example, 
p. _ —— 4i#{Bu + 2By + Br} on, 
is (Ay — Aa) (nm + 1B), Bo — Bi) ‘ 


with similar definitions for P,, , and in general for P,,_, ., . The B;, are the coefficients 


If we make use of the values of B, 
we may write the P’s in convenient form as 


in the probability density in terms of the means 


and covariances (20) and (21) 
> 42 | Ait Agly\2 9 
1 2 =F 3 | ( > s “ dt (31) 


with similar definitions for P, 


It can be shown rigorously that the characteristic function (26) has a limit as n —~ 
and that the value of this limit ‘s 
/ \ 
4 exp | ai ) Uj.l } I] V;. 32) 
\ 1 k= 
where 
~L 
| , . kra 
yi = | [A(az, t) — S(x)] sin dx (33) 
La n L } 


9 


lim I] y. = IT OP (34) 


n—« 


The details of the limiting procedure are not as obvious as they might appear on 
97 


ee | 


) 


first sight since the quantity n occurs in each term of the series and products of 
4. Derivation of averages. In this section we shall be interested in obtaining certain 
averages indicating the behavior of our system of n harmonically bound particles. Once 
we have the averaged behavior of the finite number of particles, it is but a single step 
to the case of the average properties of the string. 
We first derive the mean square deviation of the displacement of the k-th particle 
from some initial position, that is we wish 


Pi ee ee ey . = Qe 
(Vk) average — ((Z, Xk) Average T (3: ) 
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From (10) and (17) we may obtain 


ee 2 = (212) e _krr_ 
ae eS > es oul” i (36) 


With the use of (36) we obtain 


T= > (B()B,()2r0t. + cB, (1)2,oF 40 Sin w,!) 


(37) 
+ be (wre *'”?B(b)t,0P 9 Sin wt + wr twr'e "PoP .o Sin w,t Sin w t) 

where 

Big mi (cos wit a sin wt) (38) 
WW), 
and 
_ J4r(n+1)° .. re Sat 

ow A int 5g - Sh (39) 


Similarly we may derive 


= _ Ep ~ sin® kre (n + 1), Se 
wh ee »» sin” rr/(2n + 2) " Un (40) 


where 
U, = ra, (20, sin 2w,t — 8B cos 2w,t + 16a sin® as) (41) 
and 
,  2B«T r( 1)° 
E = : 3 a= < <a = (42) 
Mu ul 


Putting (37) and (40) in (35) we obtain the result (27),verace that is the mean square 
deviation for the system of n particles. 

The quantities x,, , P,) represent initial values of the displacements and momenta 
of the n particles respectively. We now average over a canonical ensemble of systems 
having all possible initial displacements and momenta after a long time has elapsed. 
The probability density P, which is a function of time, as we have found it in (21) 
approaches a limiting value as 4 > which we may call the stationary density. This 
probability density is given by 


R exp (-8 (Pio + AireTio + °° yf (43) 


Using this stationary density or as it is usually called the density corresponding to 
the canonical distribution, we may examine the following question. 

Given that a certain particle starts with the displacement @ at time zero, what is 
the mean square displacement of this particle after time /? 





ww) 
bo 
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It is seen that what we want is the evaluation 


| | exp —£ P + NiNoAIo + )pav 


The element J eo ( to be the element of volume in phase space dP i}? 
ax vee Ou ndition that the initial displacement be of magnitude 4 ts m 
estrict it l l he integ tion over the P? but does put a restriction on the integration 
I ] 
Ove | X i! ( Ol s that 
a<)> ; } é 15 
— 1 


We must thus « 


[ 4,° 
ti) 
— 
| | CXp » A \ \ 
1,’ ound 
where \ enott eleme! ( I te this C st neorem 
or I] } t) VW rie I’ ul 
i t Gy; ( 
} / G l ‘ 
) ) j 
- - - i | 
where 
I ; S s 
Ls ie 
and 
>) 2 
(,.(f — * in 1c) 
We Its fe f n particles but a rigorous 1 ical 
discussio1 process of re ng each term he limit d 
rep] ? | iB, \ at Nn tO fl le Mh nu 
ob ine ¢ nose © ned rile Her { Van Le ul ) Vl en CK 
and () r eiyv aitters ¢ od \ | i) Ol the wpre (y en 
here is t the resul 1 the finite system as well as for the string are obtained 
(ne |] rie | lon if | ec DOVE esults Cal be m ve DY recognizing that 


50) 


iad ro Sin Ara ’ 1) because U we use t I 16) 


becomes an indet t f the form 0/0. Inst 1 we therefore ust 15) and then let « > 
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so that 


: 1 “sin” krx/(n + 1) ” 
,=-5 > >4Coer'. 51 
f 2La” <4 sin’ rx/(2n + 2) (51) 


An evaluation of this yields 


p, = STB (52) 
La 


Thus the first term of (47) becomes 


E - Ek(n — k + 1) (53) 


"lie 2La°B 
and since 
. 28«T ~ 
E = His (54) 
m 
we obtain in the limit as n ~© since k/(n + 1) = 2/L the term 





‘LT P a 
UT (z _ =) (55) 


This appears also in the results of Uhlenbeck and Van Lear for the limiting case of the 


string. 
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ON THE DISTRIBUTION OF ENERGY IN NOISE- AND 
SIGNAL-MODULATED WAVES 


II. SIMULTANEOUS AMPLITUDE AND ANGLE MODULATION* 


1. Introduction. In Part I' we have determined the distribution of energy in 
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a 
is amplitude-modulated by a band of normal random noise or by a noise 
an earlier work the corresponding problem for frequency- or phase- 
es was considered. With the help of these results we now examine the 
vhich is simultaneously amplitude- and angle-modulated by a pair of 

n nome 


waves. Here, either noise wave may be delayed by an amount 


to the other, but coherence between the two modulations is assumed to 
igation of this problem is of considerable interest in the study of the 
vpical magnetron tubes, on the basis of a macroscopic model that 
netron generator as producing a high-frequency oscillation, which in turn 
slv amplitude- and phase- or frequeney-modulated by a comparatively 
nd of normal random noise, itself generated in the oscillator due to 
so-called primary noise which is inherent in the structure of the 

ing under the influence of the applied electric and magnetic fields. 
blem also provides us with results useful in communication theory when 
e appears as an amplitude and frequency distortion, or when the modu- 
n be used as a simplified model of a speech signal.t The central feature 
as distinguished from those examined previously in Part I, is the fact 
lating noise disturbances are now correlated with each other. The result 


ig” of the separate modulations, which in turn produces an asymmetrical 


im, unlike those obtained before. This is a significant phenomenon 


bsent in the simpler situations where there is no correlation. Sym- 
s can occur for selected values of the relative delay ¢, between 
nd angle-modulations. Moreover, even in extreme cases, the clipping 
lulation causes comparatively little spreading of the spectrum and 


mainly governed by the degree of angle-modulation. Special attention 


- example of coherent modulation without over-modulation. Specifically, 


rier is represented by 


9 


7 Vit (2.1) 


1951. The research reported in this paper was made possible through support 


= A,(t) exp [Zot + ¥(d)], 


(wo = 2rf,). 


19 


suboratory, Harvard University, jointly by the Navy Department (Office of Naval 


| Corps of the U.S. Army, and the U. 8. Air Force, under ONR Contract N5ori-76, 
, Quart. Appl. Math. 9, 337 (1952). 


, Technical Report No. 99, Cruft Laboratory, Harvard University (March 1, 1950). 


for the moment that any background noise accompanying the modulated carrier (and 


1 with it) is quite negligible. It is not, however, difficult to include this additional effect 
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2. The correlation function. For the amplitude modulation V,4y(é) one has from (2.7) 
of reference 1 


A,(t) = (—A,/2z) | z exp (iz[1 + kVanl(t + ts) ]) dz, (2.2a) 
JC 
and for the angle modulation one can write 
r / f° 7. 
Voar(t) = DoVen(d; Welt) = Dp | V pv (t’) dt’. (2.2b) 


As before,’ k (>0) is an amplitude-modulation index, A, is the peak amplitude of the 


unmodulated carrier, and D, is the r-m-s phase or (angular) frequenecy-deviation per 
age (or current), viz.*: D5 = (W(t)*).-av. (Ven(l vv = 


unit r-m-s modulating volt 5 = 
D: 6;/(Vrx(t)*)s-av- respectively. For the moment we distinguish 
Viav(t + ¢,) and V,p,(¢), the former of which repre- 


AanN\e TF 


wa/(V ev(t OF 
between the two (real) noise waves 
sents the disturbance producing the amplitude distortion while the latter causes the 
phase- or frequency-modulation. The effects of intervals of over-modulation, during 
which the oscillator is cut off, are accounted for in the general expression 2.22), which 
when kV4, exceeds —1 and vanishes when AV, < —1. 


gives A,[l + AVay(t + ¢,)] 
Here ¢, represents a phase lead (¢, > 0), or lag (t; < 0), of the instantaneous amplitude 
I if | 


modulation with respect to the phase- or frequency-modulation. Such a phase difference 
is included because in general the mechanism producing the coherent modulations may 
not act instantaneously: there may be a definite phase lag or lead of the one over the 
y case that we are dealing with proper modulations, i.e., those 


\ 


other. We assume in al 
whose highest significant frequency components are much less than the carrier frequency 
fo , and that a direct modulation rather than a mixing process is responsible for the 
perturbed carrier. [See section 2, Part I, for further comments on this point.| 

The auto-correlation function of the modulated carrier is written* 


R(t) = (1/2) Re{(V(to , ts) V(to + t, ts)*)s.0r. } 


= (A5/2) Re\ exp (—twot)(4°) | ze” dz | t*e"* dt (2.3) 
\ JC JC 
- [exp (tkV,z + tkV2E + OW, — TW) Jorar av. 
Here V Van(t) + t,)) is the noise amplitude-modulation at an initial time ¢, and 
V. (=Varl(t, + t& + 2) is the same disturbance at a time ¢ later; Y, and W, are re- 


and W(t, + #), the noise angle-modulation. The statistical average is 
Y, as indicated. They are not 
(,) and Vpy(t) are assumed to 


spectively V(t 
performed over the four random quantities V, 
generally independent random variables, since V4y(é + 
come from the same source and to contain in common a significant fraction of their 
total spectral distributions. We note that 


*The symbol ( ),.ar. denotes the statistical or ensemble average over the appropriate random vari- 


ables, while ( ),,. indicates the infinite time average, cf. ref. 1. 
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[exp (tk Viz + tkVoé + 1H, — (We) Jotat.av. 
= Fit, , & , ts, & 30 = characteristic function of the random 


variableskV, ,kV., WV, , and ¥, , for which 


o1 = 2, o = &, fs i, 5% = 1. 


- . e ° 3 
Since V, --- WY, are normal random variables we can write at once 


5 di :; Ine ; 
F.(¢, ’ Ce »§3,64;0 = exp (-} ys pm ttn) - F,(, é, i, a | t), 
? l 


when tf, =2, fo = &, é; = 1, i, = —l, 


where the sixteen variances u;;(j,l = 1 +--+ 4) are 


_ Mee = (ge - k°(V2)e.ev. = kv ,(0); 
My = Mai = kV, fa = kv ,(t); 


e — /yy? _ /y,?\ oe 2 ‘ 
re a (Wi) sav. _ (WP a)e.an. _ (v nn ’ 


(2.5) 
fa " fa > (—WVV2),.av. ; ais ~*~ ha kV Wi) s.09. ; 
he ta = —k(V Wo) e.av. ; Mes = Kae = k( VW). .0v. ; 
Mog = Mao = —k(V,), a 
The characteristic funetion (2.4) becomes accordingly 
F(z, &;1, —1;0 = exp [-240°'(2 + &)v.0) — ketya(d] 
(2.6) 
-exp [— DjQ(t)]-exp {—kD [za () + tA (d]}, 
in which 
A’ ’(t) A(t, t.) = (V,-(% — WVe))e.ev./Do ; 
(2.6) 
and (Vo-(V, — WVe))s.ar./Do = A(t) = —A(—t#, &), 
and Q(t) is given by’ 
Q(t) = [ Waal fC — cos wt) df, or (2.7a) 
= [| Wr — coset) aff’ (2.7b) 


Jo 
for phase- or frequency-modulation respectively; ¥4(t) is the auto-correlation function 
of the noise wave V,,(¢) which produces the amplitude-modulation. 

There remains now to evaluate the cross-correlation function A(t, t,), relating the 
two modulating waves. We consider first phase-modulation. As in previous cases we 
express the modulating waves in terms of their Fourier transforms S;(f),4 and S7(f)r 


3—D. Middleton, Quart. Appl. Math. 5, 445 (1948). 
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for the amplitude and angle modulations respectively. Replacing statistical averages by 
their equivalent time averages, in virtue of the assumed stationarity (and ergodicity) 
of the random processes, we have finally 


A(t, t,) | lim 7 '{S, f)aSr(— f)r} exp (ctw)(1 — e**’) df. 


n terms of a modulus | S,(f)4 


Writing the amplitude spectral densities S7(/)4 , ete 
(= | S7(—/JS)4 |) and a phase ¢(/) 4 o with similar expressions for S;(/)- , 


since we are dealing with a real wave, and noting furthermore that the mean spectral 
intensities of V4,(¢) and | f) are defined by 
i W pv(f) lim 2|Sr(f)r |T , 


7 


we can write the cross-correlation function A as 


(2.9 


t+ bd, — br) — COS lwl(h, — bt) + ba or]} df, 


~ ' 


nsities are even functions of frequency, while the phase factors Q; 


since the spectral densi 
by 


and ¢, are odd-functions. In the case of frequency-modulation we replace S 
S,(f’) iw’ in the above, according to Eq. (2.2b), and obtain finally the cross-correlation 


function 


f ‘ 
‘7 i | f 4 P\] 
Arm Ee | } 1] PM } | 
, 2.10 
fain (wl, +- o, — dr) — sin (w(t, — 1) + b, — or)} Af /w. 
The cross-correlation function A(¢, ¢,) defined above in (2.6a) has a number of im- 


portant properties 


7) We observe at once from (2.9) or (2.10) that if the two noise waves have no 
spectral region in common, 1.e., W4, und Wp y(f) do not overlap, the cross-correlation 
function vanishes, even though the two disturbances arise from a common source Thus 
the amplitude- and frequency-modulations are no longer in any way correlated. In 


have a significant proportion of their power 


general, however, Vi and Wyy(/) may hav 
distributed in the same spectral region; the precise amount is measured by a quantity 


which we sl il] eall the 0 erlap- Ol cohe en funel on x c.). V1Z.° 


x(t) = (VW). or. /Do(= (V¥e)s.0,/D 
| [Waun(f)W f)| cos (wl, + bd4 — or) df 


(2.11a) 


—] [Wax(f)Wen()]'” sin (ot, + ¢4 — or) df/w 


(2.11b) 
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respectively for phase- or frequency-modulation; A, as we shall see presently, cf. (2.15), 
is a constant depending on the behavior of both spectral densities W4y(f) and Wery(f) 
as 7 — Q. 

(77) When ¢ approaches zero, A(¢, (,) vanishes identically; this is a direct consequence 
of the definition (2.6a) of the cross-correlation function. 

(777) When the time lag (or lead) ¢, of the amplitude-modulation with respect to 
the angle-modulation becomes indefinitely great (i.e., when t, — + ©) we observe in 
the case of phase-modulation that since 


lim [cos (wt, + ¢d4 — dr) — cos (wl, + b4 — br — wl)] 


t 


(2.12) 
= F2zw sin ef 5 — 0) = +e sin at). 5p - @, 
we have* 
‘ . l , , 1/2 . wl 
lim Ag wd: t,) = + lim ‘3 [WW an( fi on(f)] ® Sin 9 \ (2.13) 


which usually vanishes, unless the spectral densities together approach zero as O0(w "), 
n > 2). For a band of noise whose zero-frequency density is finite or zero, (2.13) always 
vanishes. In a similar fashion one finds that in frequency-modulation 

. l a . “ a 2 - 

lim Apw(t, 4) = 4 | (Wav(NWev( fl [6(f — 0) — 6(f — 0)) df = 0 (2.14) 


2 J. 


identically and quite independently of the spectral distributions of the modulating 
noise waves in the neighborhood of f = 0. The vanishing of the cross-correlation function 
in the limit of infinite delay times is explained if we note that the original wave-forms 
are now so separated in time that even a partial coherence between the two waves is 
completely lost. This is not true however when the major portion of the noise waves’ 
energy lies in the almost zero-frequency region, ef. (2.13), as then the waveform changes 
so slowly in the course of time that even an infinite delay does not destroy coherence 
entirely. This can also be seen directly from the overlap function x(/,) when (, —~ +. 
iv) A similar argument may be applied when | ¢ | — © to show that 


* « 


t 


lim A wif, t Koa( ty) 3 lim Apalt, t,) = Xra(ts) a ; Wv(0)W py(0); (2.15) 


here the time-dependent part of the cross-correlation function may or may not vanish, 
depending again on how much of the wave’s energy is located in the spectral region 


U, e), as e — U. 
v) If the amplitude- and angle-modulations have the same spectra, i.e., Wyay(f) = 
Wrev(S)3 ba é@p , a simplifying assumption that is ordinarily not at all critical, the 


calculation of the auto-correlation function (2.3) and the corresponding mean intensity 
spectrum W(/) is greatly reduced, since the cross-correaltion function A takes a com- 
paratively simple form. Note also that when the delay ¢, vanishes, one has 


Agu(t, 0) = Qsw(d), cf. (2.7a); Wall) = Qsy(o) — Q(t) (2.16a) 
*A rigorous treatment recognizes these results as a form of the Dirichlet conditions.‘ 
4H. S. Carslaw, Theory of Fourier’s Series and Integrals, 3rd rev. ed. (Dover, N. Y., 1949), pp. 92-94. 
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Aru(t, 0) = Qru(t), ef. (2.7b);2 = Pal) = Qra(i). (2.16b) 
For example, consider the case of identical gaussian noise spectra, W(f) = W. exp 
[—w’/ws]; one readily finds from (2.9)-(2.11) that 
Agw(t, t) W ow, /42'’*) [exp (—wits/4) — exp (—w;(t — t,)°/4)]; (2.17a) 
Xou(ts) = (Wow,/4r'’”) exp (—wyt/4), 
Aru(t, ts) = (Wo/4)[A@t./2) — Mw@lt, — t]/2)]; 
(2.17b) 
+ 9 sel te) W./4)A(w,t,/2), 
where @)(x) = 2r°'”” [5 exp (—y’) dy. 


The complete auto-correlation function (2.3) for the noise-modulated carrier can 


now be written 
R(t) = (Ao/2) exp {— Dp 2(0)} Re exp (—is t)(47r°) 


| 27? dz exp [iz — k’az?/2 — kDA“ (Oz) (2.18a) 


dé exp [7% — kD. A (HE — kt /2]-exp (-Kv.(eI} 


ee 
" 


= (A3/2) exp {—D2(t)} >> [(—k’va)"/n!]ro(0) 


» YY (H—kD (dade) Oy OWA (! (2.18b) 


» Refexp (—twot)ho.nsjRo.nsi}s 


where 


ho.esn = (2x) | exp (iz — kK’ yz" /2)2"*"~ dz, 
J © 
and (2.18c) 
A (b)/(babr) =A (OH; Wa.r(O)/Pa.r = roba.r - 
The amplitude functions h are given explicitly by Eqs. (2.11-2.12) of reference 1. 


The spectrum is found in the usual way by taking the Fourier transform of F(t). 


3. Remarks on power and spectra. As before, the mean total power W, is obtained if 
above, namely in the auto-correlation function R(t), and the 
2.18b) also, as ¢ is allowed to become infinite. In 
> 0,1 = j = 0), while in the latter 


one sets ¢ = 0 in (2.18b 
mean carrier power W,, follows from 
the former case one gets a single infinite series (n - 
a double series development (n = 0, j, 1 > 0), is the result. W, and W,, may, however, 
be obtained in closed form, if we follow the procedure of section 2, ref. 1. We need now 


the joint first-order probability density W,(y, x) for the two correlated, normal random 
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variables y = kVay(t + ty) and x = Di Vyy(t) or Dy J‘ Vew(t’) dt’. This may be obtained 
directly from the characteristic function (2.4) and (2.5) for the second-order distribution 
on setting ¢ ¢, = 0, which gives the characteristic function F,(¢, , {). 


By inversion one finds® that 


, l = ; ;' ta i 
Wy, x) = 3 I| Fi(gi , &s) exp (—ty$i — txfs) df, dfs 


(27r) 


bo 


rA'’*)~* exp {34 [DjQ(@)y? + kar? — 2Daxitxul}, (3.1) 


A = k’?W,4D5Q(e) — (kDox)’. 
Expressing (4), (2.1), now in terms of the random variables y and x, namely, 


V(t) = Ao(1 + y) exp (2x + tal), (y => —1); = 90, yy S$ —}), 


W, = RO) = ; {| A,(te"*” Pica 
= (43/2): [  dyaty? | "Wily, 2) dx (3.2) 
vy=-!1 -o 


= (A,/2) (i + y) Wi(y) dy, 
a | 
where HW’,(y) is a normal probability density for which y = 0, y° = 4k’, showing that 
the mean total power is the same as if there were amplitude-modulation only. This is 
what we would expect, since angle-modulation merely redistributes the original wave’s 
energy among the carrier and sidebands without changing the total power (see sec. 2 
2). On the other hand, the mean power W’,, remaining in the carrier is influenced 


of ref. 2). 
by the phase- or frequency-modulation, because of the coherence between modulations. 


We have 
, l wV(t)y 12 2 /é | . hi rT |? 
Wy. = 5 | (Aolde™ sue. ? = (43/2)| | A+ydy] e*Wiy, 2) dz | 
“ a | J —@ 
= (43/2)| [ ata | en hd, Dae, 
(3.3a) 
™ ] 1 a @ 
= (A?/2) exp [— D32()] | =i | dz 
| 7 J —(2k9ya)—3/8 


12 


“(1 + 2[2k*y4]'*) exp [—@ — 1Dox/(2¥4)'")"] |, — (3.3b) 


which may be integrated in straightforward fashion with the help of the MacLaurin 
series for ¢°' (2 + a). We find the mean carrier power to be finally 
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W,, = (A?/2) exp [— D52(@)] | ho.o + tkDox | hor | 


2 


1/2 ~ l ° 1/277 n-2) 2 1/2 
+ kW? Do eDox/ Wa Vo" (val) |, 


n 


representative values of which are shown in Fig. 3.1, along with curves for the mean 
continuum power IW W.-W Typical curves for W, are given in Fig. 2.2 of ref. 1.) 
It is immediately apparent from (3.3b) that if there is no angle- or amplitude- 
modulation (i.e., D, — 0, or k — 0) one obtains the previous result for W,, , namely 


Eqs. (2.21), ref. 2 or (2.16), ref. 1, with corresponding expressions for the mean con- 


tinuum power WW’. . Furthermore, if there is no coherence between the modulations (so 


that x — 0), the resulting power is proportional to the product of the individual carrier 


yowers for angle- and amplitude-modulation alone, viz.: 
I | 


W 15/2)he,o-exp [— DoQ(~)], x — 0). 3.5 


We distinguish now between phase- and frequency-modulation, according to the 


remarks of sec. 2, ref. 2. Figure 3.1 shows only the case of phase-modulation by noise 
with a gaussian spectrum [ol which W.,,(0) is finite,* ef. Kq 2.35), ref. 2. Increasing 


6 (= EV (h)« ), either by lessening the rate of sweep or by increasing the r-m-s 


phase deviation D, , decreases the amount of power remaining in the carrier and makes 


noticeably greater the amount of energy available in the continuum. However, because 


le of this effect is determined bv 


of coupling with the amplitude-modulation, the sca 
a), the extent of the amplitude variation (—y,A"), and (b), the amount of coherence 
(~ x), which depends jointly on the time-delay /, and on the magnitude of the angle- 
modulation (—D5y For very heavy over-modulation Eq. (3.4) reduces to 
W,, = (Ao/2) | [ko 4 /2r 1 + exp (Dix /2W4) — (2° Dox/Wa- 

| + 4 - 1D xT 4 4 i 

exp j ay + , 172 o.U 
(24) 
exp Di, Q()], k—o); 


The integral in (3.6) is tabulated on page 32 of Jahnke and Emde’s Tables (Dover, 


1943).° Integrating 3.31 DY parts o1Ves the appropriate asymptotic development when 
Dox/(2W4)'“~ is large, for any value of k(>0). We have finallyj 
W,. = (A5/2)(k’ Wi /2rDox’)-| (1 + 2¢/kDox + +--+) 
9-7 
oO. 
exp {—1/k°W4 + Dilx’/va — 2(@)]}, (Dox/(2Wa)'”” > 1). 


We return to the expression (2.18b) for the auto-correlation function of the modu- 


lated carrier. Applying the theorem of Wiener and Khintchine [(1.6), reference 1], we 





*Since Wry(0) > O for the same type of spectrum, Qry(@) becomes infinite in such a way that 
W;, above vanishes. Consequently, all the wave’s energy is distributed in the continuum; therefore 
W, = W. and W, are independent of x and Dy. 


5Jahnke and Emde, Tables of Functions, 4th rev. ed. (Dover, N. Y., 1945). 
tk never appears to any lower power than k~ in the series in Eq. (3.7). 
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distinguish three general types of spectral contribution: (i) ‘“‘pure’’ angle-modulation of 
the carrier by noise, i.e., (noise X carrier), or (n X c) terms, for which (n = 7 = 1 = 0); 
(ii) “pure” angle-modulation times ‘‘pure’’ amplitude-modulation harmonics, namely 
(n X Cencie’ (2 X C)am components, for which in turn (n > 1,7 = 1 = O); and finally 
(iii), the cross-modulation products (n X ¢)asie KX (" X Caw, (n 20,7 +12] JD, 
which arise from the coherence of the simultaneous angle- and amplitude-modulations. 
(We use the term “pure” above to indicate the independent or incoherent character of 
the modulation vis-a-vis the carrier and the other modulation.) Components of types 


(i) and (ii) are readily seen to yield only the usual symmetrical power spectra (about 
FV Modulation products belonging to group 


yield both symmetrical and anti-symmetrical distributions. As we shall see presently, 


a direct consequence of the coherence between the two noise 


ill), however, are more complex: these 


this anti-symmetry is 


modulations and depends critically on the amount of delay ¢, between them. 


We observe that the terms in (2.18b) for which 7 + 1 odd (n 1) are responsible 
for the spectral asymmetry. To see this we note first that terms involving A and 
\‘~’ (2) in the auto-correlation function R(t), that are odd in ¢, yield spectral contributions 
which are odd (i.e., antisymmetrical) in . Furthermore, such terms have as a 
factor sin wf, so that R(/) itself remains an even function, as required from the nature 
of R(t). Now ) and A may be split into two parts, one of which is even and the 
other odd, in such vy that if we let u [) l j!l! represent the summand of 


the double series in (2.18b), including the amplitude factors hy, , we may write sym- 


bolically E, +0 EE, + 0 i!, where # and O are respectively the even and 
u(t). Then by direct expansion one gets for + | odd) only odd 


odd parts 
terms, and fo +- ] even), only even terms. Let us illustrate with the simplest, 
but most important for which rd ks 0), and to facilitate the discussion 
without i1 ny way changing its essential features, let us further require that D y be 


so small that we can set exp DoQ(t)] equ: 


to unity We have as before to distinguish 


between the generally dissimilar situations of phase- or frequency-modulation. From 
2.18b), (2.9) and (2.10 e have for the part of the auto-correlation funetion for which 
§, +. | 1.7? Q 
R(t), = (A5/2)kDo(Waver) [A (8) + ASO] Refexp (—twot)ho.o-ho.1} 

3.38) 


6M) 
FMs 


, 


sin (wt, + 64 — op) 


‘ df’ 3. Sb) 
COS (w’ ft — oo ~~ op) | 
respectively for phase- or frequency-modulation. The power spectrum for these cross- 


products becomes 


Wf) am.rm — AskD i} ye 


+1 Qemem(S 16 — [f+ f) - «f -[-f-f) (3-9) 
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where 
Qsulf’) = (Wav) Worf)” sin [w’ts + Of — oF], 

and (3.9a) 
Qeulf’) = (Wav) Wen(f)]'” cos @'t, + 64 — o7)/o’. 


Since both /, and f are positive, for obvious physical reasons, it follows immediately 
that the spectral density (3.9) reduces to 


Af—-f), hhr>f>o 
W(f) = A2kD,(—ilo.o-Ro.1) (3.10) 
lay —f.), f>&h>O0 


inasmuch as the integrals over Q in (3.9) are zero for all negative values of f, and f, 


and because ( 3) is \ anishingly small, Q) being a low-fre quency spectral distribution 
of the type W, exp w /w,). The asymmetrical nature of the spectrum is at once 
evident from (3.10). Actually, the discontinuity at f = f, does not exist, but is due 


here to our simplifying approximation that exp [—DjQ()] = 1. Including the expo- 


wif | 9 =CONSTANT Ww (f) W, /W, = CONSTANT 
d (FMAM) co} b 











a + —— - 
Y 
c hs W,!, 
p! { AM LEADS ANGLE - MODULATION 





) (Ww,t,>0) 
1 Sol , 
(a) 


Fic. 3.2. Symmetrical and asymmetrical noise spectra for correlated noise modulation. 
nential factor removes the discontinuity but does not alter the asymmetry; a specific 
example is considered in section 4(b) following. Similar remarks apply to the terms in 
R(t), (2.18b), for which n >1 and 1 + 7 = odd >1. 

Figure (3.2) shows typical spectra for a fixed amount of over-modulation and for 
various values of delay; (a) illustrates the case of phase-modulation by noise, while (b) 
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| 


gives the corresponding case of frequency-modulation by the same noise wave. The 


noise causing the angle-modulation is here identical with that producing the amplitude- 
distortion of the carrier, but is delayed by /, seconds. We note at once that phase- 
modulation /eads the frequency-modulation by approximately a quarter period. This is 
understandable if we remember that the frequency-modulation is the infegral of the 


ph: se-modulation (for the same modulating wave), cf. (3.8), and consequently lags 


behind the latter in time by a delay t, — m/ 2w, on the average. In the case of simultaneous 
} 


amplitude- and angle n odulation by i sine-wave alone this 1s exactly the delay Now, 
] 
j 


as the de lay or lead becomes progressively greater, the maxima of the spectral distribu- 
distributions themselves, oscillate about f, in ever-decreasing 


tions, and therefore the d 
} 


swings ul for infinite lags or leads all coherence between the modulations is lost, and 
only a symmetrical spectrum remains. The decay of the swings’ maxima is quite rapid 
for phase-modulation 1 somewhat less so for frequency-modulation, ¢f. (2.17 rhe 
oscillations of the spe maxima are also shown,in Fig. (3.2) above; we remark, how- 
ever, that values (of at the points of symmetry (/ f,) and maximum asvmmetry 
are only approxi! nd do not represent a purely periodic phenomenon. The precise 
values depet i in COI licated Way on the spectral shape of the modulating noise, the 
amount of over-mor tion, and the intensity of the angle-modulation. From this we 
see also that a symme cal spectrum alone is not sufficient evidence upon hich to 
postulate lack of e nee between modulations. For a given level of clipping, the 
power (1.e. the re nae the spectral curve in the continuum is a const: nde- 
pendent of the phase ferences between modulations, (but not necessarily independent 
ol the type of angle-modulation 
As noted above, the en rey spectrum may be resolved into three princip: | orders of 

modulation products (i)-(ii1), viz 
Wf) = WC + Wf 

>. 11) 

i Wf), Wifi + WO 

respecti\ When the deviation of the angle-modulation is very large compared to the 
rate of deviation, so that asymptotic developments are. possible in the manner of 
section 2 (reference 2), we find from (2.18) that 
Wren AZ 2, va)" /n'}) 

3.12) 

bh | c* (Oro(t), exp — Di 2’ /2) cos (w, — w)t dt 

where the polynomials « t) yield the asymptotic development in descending powers 
of w,/w, or 6; ; thus t) yields (wz/a,) ~* or 6;°* multiplied by appropriate constants 


Specifically we have 


; me me . ' -~Doo” ., Ba” 
ss weg » = ( oY omr fF 


—D2 , 2 oe 
(6) a a | | 4 4 10 2 23. 
ev = ( gy ° 216141 2% -* gigp 3 
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wa {— Dia” no Di | a" it ol 
. ( gs (- 10! =< o 2! 6 + - gia! t 
(3.13) 
|. ie i Dia* ., 
gure. ¢ t+ ara’ ji ete 
where 


Qe im (~ 9)" G0 ; n > 1; OM = b.f-1)""; n > 0; 


ox 


b, a w W(f)cea.em df. 


The spectrum of the ecross-modulation is more complicated, but can be handled in a 


similar fashion. One has then to consider terms of the following type: 
} F n ) i / l 
-1) Re) .| | ri (t)"e?? (DAL, t)' ACE, te) 
- exp (— DiQ"”’t/2 — iwot) cos wt at} 


he amplitude-modulation is excessive, so that over-modulation 


On the other hand, if the 
curs approximately half the time, the asymmetrical components are suppressed, and 
one has essentially a symmetrical spectrum, which nevertheless is influenced by the 
coheret een the angle- and amplitude-modulations. The analytical explanation 
imme nly / n + q = even) does not vanish when (2/°y,)'~ >, except 
/ , whi mes independent of (24°y,)'*, ef. Eqs. (2.12a-e), ref. 1. Consequently, 
is alv even, and only symmetrical terms remain. Since these components are 
| ( L,), the contribution of the single remaining antisymmetrical pair 


0), which are O(kW'/"), becomes ignorable in the limit (24°y,)'* ~@, 
the amplitude distortion which effectively determines the “cut-off” 


oscillator. The angle-modulation alone does not alter the amplitude, 


spacing between zero-crossings, and therefore one expects that the 
‘on’’- and “off’-periods of 


modulation controls the duration of the ‘ 


trong amptiil le 


the wave 

4. Coherent and incoherent modulations. ‘wo important cases remain to be consid- 
ered in some detail: one, the situation of no coherence (A = 0) between modulations, 
and fwo, the case of coherence, but no over-modulation. A discussion of one and two 
follows belo In parts a) and (b). 


The general problem simplifies greatly when there is no coherence or coupling be- 
tween the modulating waves, a condition represented analytically by the fact that 
anishes, and physically, by the complete independence of the noise waves. The 


A , l 
auto-correlation function (2.18) reduees now to 
R(t) = 7,(t Ro(b ay COS Wol 


(4.1) 


= (A2/2) exp [—D22(t)] SS {((—ka)"/n!ro(OAh5.n COS wol, 
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=W(f); + Wf), ef. Eq. (3.11), contains no cross 


and the intensity spectrum JW (/) 
terms of the type W(/ . The mean total power in the wave is given by Eq. (2.17), 
reference 1, as before; the power in the carrier and continuum is also easily found from 
(3.4). 
Let us consider the specific example of phase-modulation, and let us assume a gaussian 
spectral distribution for the modulating noise; (in general, the two noise waves need 
not have identical spectra or powers.) The spectrum follows at once from the cosine 
Fourier transform of E 1.1) and we obtain finally with the help of (c), section 2, 


reference 2, 


ar ior” exp 6 ; oe 6°" exp (— 8°/¢q) 
VW (}) = : ) No oof — Jo) “Tr # bs P 172 u 
w \ q=1 q-¢ 
(4.2) 
Wy f J , ( 3 y - \* > ) ) —_ 

. he hwy exp D /\G TT MWs4/e 5 _ @®o wd 

jw 2 ' ' , 2 2\1/2 € p = 
: Q ? (q + nw4/w>) J ws, 


The first term represents the carrier’s contribution, the second, the spectrum qaue to 
attributable to 


phase-modulation alone, and the final term, the distribution W(f),, 
the noncoherent modulation products between the angle- and amplitude-modulations. 
When excessive over-modulation occurs, we find that 


; ( 2 2 
= - Ws. exp 0 : " . = 5 4 0 TR Mag 7) 
Wifes y Lf- f+ & (: p(- 8/9 


m exp |—B'/(q + wa/ws)] 
2(¢ + w/w) ) 
O(n) exp [6° + (n + 2)w? /a)] 
“ \n!?2"(2n + 2)(2n + 1)[q + (n + 2)w2/wr]'” 
t.3) 
2n)! exp (—8°/(n + 2)w4/w;) ))! . 
T nl’2 2n + 2)(2n L)[(n + 2)w4 /w | })f» _. 
and when there is ignorable over-modulation we have simply 
{ 
exp ff | A- 
= 5 T exp (—p /q) 


W(f) = Aj” if —- f)+ LD 
2 q°9 


— — @ exp |—f gq+-w "| Ww, ° 1/2 : 
+ K Wa rs ; D172 ; 2k Wa) < 0.6. 
qf: / “T- Wa/W,/) 
Figure (4.1) shows the energy spectrum for (2k°y)'’~ = 0.5 and for various values of 


very slight broadening of the spectrum results 


ratio 6 \ general, but 
de- 


the deviation 

from the additional modulation. Similar remarks apply for frequency-modulation 

tailed calculations for which may be carried out along the lines of section 2(b), reference 

2). In general, because the spread in the spectrum due to clipping is small, even in the 

extreme case of 50 per cent over-modulation, as can be seen from Fig. (2.1), reference 

1, we expect little additional broadening of the spectrum due to this effect. For large 
' is obtained in asymptotic form from 


values of 6; or w;/w, the spectral intensity W(/) 
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(3.12) and (3.13); in the instance of a gaussian spectrum for the amplitude-modulation 


(AW, exp (—w/w,)) we have 
—. (—k’ Wa)" 22 ’ ' oe 
wne> ae h2_, {Eq. (2.33), ref. 2, modified}, (4.5) 
where now Kq. (2.33), reference 2 is easily modified by first multiplying numerator and 


denominator of each term which contains a factor ¢* by a suitable power of Dj , so 






































“ T T T T l l l T 
—— — A FR 
I = W(f). = MEAN SPECTRAL INTENSITY wif), + W(f)p “7 
1 FOR PHASE -MODULATION ALONE : 
8 
t 
' 
a 
mI 
1 | | | 
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| - 
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Fic. 4.1. Mean intensity spectrum for simultaneous uncorrelated amplitude and phase-modulation of a 
carrier by noise (no over-modulation). 


that [D52"°’]’ results in the denominator, and then replacing [Dj2°*’]’ so obtained by 


[Dian + win/2]*; the argument of the ¢*” is similarly changed from 8” to 


m\ ys 


B’'/(1 + win/2D,2”)™”. 
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(b) Coherence, but no over-modulation 

This is a particularly important case because it is the simplest one which incorporates 
all the new features of coherence between modulations, and because it offers a promising 
macroscopic model for the behavior of noise in a magnetron, as described at the be- 


ginning of this section 
When over-modulation effects are ignorable (<1 per cent of the time), only the 


amplitude functions / =1) and hy, (=2) yield significant contributions, and we 
obtain from Eq. (2.18b) the following spectral components: (noise X Carrier) gacie mod 

for which (n = 7 = 1 0); (n XK Chante’ (% X Claw , for which (n = 1,7 =l1 0); 
and finally (n X « xX (n & Chay, Where (n = 0,1 + 7 = 1), all of which are repre- 


sented in (4.6) below in the order listed: 


+ 2) D,\ | Qf’) sin w’t aj’) sin wot ?, (2h?y,)'? < 0.6. 
\ / ) 

Here Q(/’) is given by Eq. (3.9a). We observe at once that the mean total power is the 
same as if there were no angle-modulation at all, which is not surprising, as angle- 
modulati mn alone doe not « hange the energy content of a modulated earrier. Note also 
that the ross-correlation term never contributes to the total power by definition, ef. 
(2.6a)), but represent edistribution of the existing energy, determined by the carrier 
strength and degree of amplitude-modulation. The mean power in the carrier is the 
coefficient of cos r sin wf as ¢ becomes infinite: 

HI 1,,/2) exp D2 + (kD,/2) lim [w’Q(f’)}}, 1.7 
showing that IV,, is determined essentially by the type of angle-modulation. The second 
term of (4.7) vanishes for the usual types of phase-modulation, while for frequency- 
modulation the exponential factor vanishes, as explained earlier in section 2, reference 
2, with the result that H is here described by the previous expression 2.21), reference 


2. The continuum powe 
i A?/2)(1 + k’v, — exp [—D20(@)]). (4.8) 


Let us consider specifically the example of phase-modulation when the modulating 


i 


] 


ves ire identie so that 


noise w 
Wav(f = W ; f W exp (— @ Wy ; Ws = VY: W yw hor 
(4.9) 


‘Os = OF: r,(t) = exp [—w,°/4]. 
To obtain the desired spectrum we need the following integrals: 


a J sin at 
\ -exp(—Ot t ct) dt 


Leos at) 


(4.10) 


{ m > \ 
eaters — ; J sin (t+ac/2b°) | 
= (r b) exp [( — a’)/4b°] ‘ - Re(ib!) > 0. 


‘Leos (+ac/2b°)) 
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Fic. 4.2. Cross-spectral intensity of a carrier simultaneously amplitude and phase modulated by noise 
(no over-modulation). 








Ion 
Om 


i. 


ae) 


DAVID MIDDLETON 


[Vol. X, No. 1 














at 

°o 

win) oP : ; 
rs o @ ro) rn 

oa am co ae coe nen ee ee 














_ oo 
/ £ - e 5) 
y, = _ x 
a Ss rn] } 
" ~ €I 
7 VZ20 1 
/ ‘ Oot ° 
j ba - > 
/ y 7 > = Th) Oo 
f// 4 z - m w 
y - Pb 
Nfl f Oo 
wit mom 
Mod a or 
DS WH] ~~ . 
yy 
} id Db = 
| / We; = a 
| pe - = 
| fi m m 
J m 
4 
F; 
n WA » | | ! | 








Mean spec tra 


intensity 


f a carrier simultaneously amplitude and phase modulated bv corre- 


lated random noise (no over-modulation, 64 constant, various delays). 


1952] NOISE- AND SIGNAL-MODULATED WAVES 53 


az a® 
wa & wif), +W(f) + wif — 
W(f ) OQ I r ly WO, 










































' ° a) b fon) @ ° Oe) Db 
” n l T T T T T 
\ \ 
‘ \ 
\ \ 
\ \ 
hy \ 
f ‘ \ \ 
re‘ % 
- . > \ = 
\ \ \ 
| \ \ 
\ \ oun @® 
{ \ \ 
\\ \ *% 
\ \ a oe aig 
Db ' \ oO Pps + 
| 2202} 
oo os 
COO-u| OF 
\ OVBvS/F GB PC 
| | we Nizdave N 
“— se “- 
| ee = 
| <— ° 
a as 2 — 
t 
\ 
b \ _ 
| % 
j > 
| oO 
| db 
D 
o FD +4 
uo m 
/ D 
| } 
| 
] 
ce 
5 | 
} VA | ~~] 
1 | 
| 
l 
| is -# : 
Wf / > = oc 4 
| iW // "| a 
j oo +” 
| / . 
f a oe 
/ " = 
| z=wv °o en Uw 
x 
/ o> - - 
Y ev > = 
rae, be = 
| / am . 
€; of z 
= S 
. i n < 
fe7 ? oa ” 
b / r m s oa 
4 
f = = 
| / m 
| m 
'/ / > 
I | / 
, 
| { 
wn | ft > 
1] / 
yy / 
fy 
|i 
‘] /i 
f, 
/ 
6 Ul £ / | N l | i i 








Fria. 4.4. Mean spectral intensity of a carrier simultaneously amplitude and phase modulated by correlated 
random noise (no over-modulation, t, constant, various degrees of modulation). 
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Application of the Wiener-Khintchine theorem yields 


Wf) = Wfloa + Wfom-aa + WAeowxam 5 1.11) 
where JI y = VW is given by Eq. (2.48), reference 2, W(f)ou-a0 = Wf) is 
obtained dire tly from the last term of (4.4) when ws = Ws, and 
W(fleaxase = Wf) = Ok ¥ “= exp (— 86 (Aor ’* /e) 

A f { @ t.B ) 
= 7 W72 &X] °/(n + 1)\2sin (——" ) exp [—nuts)’/4@ + D]f, (4-12) 
rn: | /l Zz } 


j= w W)/W 
Typical spectra are shown in Figs. (4.2)-(4.4) for a variety of time-delays and for various 
values of 6, . The spectra are in general asymmetrical, in the manner of Fig. (3.2); the 
spectral spread increases with the larger values of @,. We remark that, unlike the case 


of no coherence, it is not possible here to speak of a phase-modulating noise equivalent 


to a frequency-modulating wave (in the sense of Eq. 


2.19), reference 2) which produces 
the same final spectrum, even if we assume the same disturbance causes both amplitude- 
and angle-modulations. The explanation lies in the entirely different types of modulation, 
which are coupled together in such a way that altering one does not produce a com- 

the other. When the angle-modulation is sufficiently intense 1.e., 


pensating change 
the needed asymptotic 


high modulation indices), we obtain from (3.12) and (3.13) 


representation of the spectrum, namely 


W(f = 4 : | eet) teJ Wa) (t),) COS (wy — wt t+ 2k-D, sin (wa, — w)t 


Sent 
t.13) 
| Q(f’, t,) sin w’t df fe (0) exp (— net 2). 
The precis¢ forms of Q etc. that appe: in the exponent and in the polynomials aie (3 
Eq. 13). are described in Eq. (2.34), reference 2. W and W'(f);, follow from sec. 2, 


reference 2, and the above W’(f),;;; has the following general form 


araee jee We eee ie 
Wie -AD A a 1.14) 


2m + 1)![DoQ” |" 


where 


and G 3) is the (2m + 1)st derivative with respect to 3 of the expression in the 
parenth ses Of Eq 2.33), reference 
with a gaussian spectrum one easily shows that when the amplitude- and angle-modula- 


2. For phase- or frequency-modulation by noise 


tions are identical 


Wf) ar3 = (27) 2k y) +s '% 
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or 
— ye As 2 2 wa (—1)” (wr \" 9°” (wots /2'"”) 
W(f) = —x(2m)'”? — [2(2k’y)]'” apoueal taal Mite gr ter.* fe 
oi Wa [2(2k'y)] p> 2” \wa (2m + 1)! 


(4.16) 


_ 


- Ge"); (w,/wa)? K p’ = (= —#)(#). 

Wp Wa 
We observe that if the amplitude-modulation is particularly weak, the (angle X A), 
or W(/)i:; , components are strong compared to the incoherent terms W(/);, , since the 
latter are O(k*~,), whereas the former appear as O(/°y,,)'. The “pure” angle-modulation 
terms ]¥V(/), , however, remains unaffected. 

5. Remarks on signal and noise modulations. When a signal accompanies the noise 
there are a number of interesting possibilities: (i) amplitude-modulation by signal and 
noise and angle-modulation by noise, (ii), amplitude-modulation by noise and angle- 
modulation by a signal and noise, and (iii), angle- and amplitude-modulations by signal 


and noise; in each case the noise waves and the signal waves are respectively coherent, 


but there is no coupling between signal and noise. We shall outline the principal steps 
in the process of obtaining the spectrum, but without giving a detailed treatment, in- 
asmuch as the more complex problems of signal and noise in this case introduce no new 
concepts or techniques not already considered in the preceding sections. 

Since (iii) contains (1) and (ii) as special cases, let us examine the more general case 


first. The modulated carrier 1s now 


qty V(t) = A(t s.y exp [tot + VWOv + ¥()s)| (5.1) 
vhere A,(! is an obvious generalization of (2.2a), viz." 
Ld) 1, /2z) | z ° dz exp fiz[l + kVay(t + 4) + wVs(t + 6)]} (5.2) 
JC 
and V(/ , wii are specified by | 24 a (2.2). 3.1b), ref. 21, or by (Zr The auto- 
correlation function (2.3) becomes 
( . 
R l A; 2 ReSe tr’) | z _ dz 
\ Cc 
| &-%e'® dt [exp {tkVmz + tkhVvak + Wy — Wyo} Jocatcav. (5.3) 
¢ 


exp {tuVsi2 + tV go + 1Vs1 — i s2} Iotat.ar-f 


where the subscripts 1 and 2 refer to the wave at an initial and final time, ¢ seconds apart. 
The first ensemble average has already been found, namely Eq. (2.6), and the second, 
which is the characteristic function for the joint signal modulations may be expressed as* 

*It is assumed that. the fundamental frequencies of the two signals are here commensurable; otherwise 
there will be no coherence between signals, and the characteristic function then factors into the product 
of two characteristic functions F(z, &; t)s.am°Fo(z, & Os.angie wv » Tepresenting the (now) independent 


amplitude- and angle-modulations. 
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=t 1 exp [iuVs(to + tee + iuVs(to + ty + DE + W(to)s — W(t + ts] dt (5.4) 


= i & e,(—1)”C,,(2)C,,,(E) COS Maral, 


since F, is an even function of ¢t. The auto-correlation function (2.18b) may still be 
used, provided that according to 
(3.7) of ref. 1, with C,,(z) in place of B,,(z) therein, insert 
tor, and perform the summation over all values of m. The precise 
of course, on the wave shape of the modulating signals, and 


very complicated functions. In general, the 


we replace the amplitude functions ho.ns+¢ DY An.n+q 
—1)"e,, COS mw,f as an addi- 


tional trigonometric fa 
form of C,,(z) depend 
except for the simple sinusoidal cases, are 
spectrum will be asymmetrical in the discrete or signal portions of the distribution, as 
well as in the continuum. 
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GRAVITY WAVES OF FINITE AMPLITUDE 
III. Steady, Symmetrical, Periodic Waves in a Channel of Finite Depth* 


BY 
T. V. DAVIES 


King’s College, Strand, London 


Introduction. Part I of this series** is concerned with the theory of finite amplitude 
gravity waves in a channel of infinite depth, Part II deals with the problem of the solitary 
wave. In the present paper, the finite amplitude wave problem is considered for a channel 
of finite depth and it will be shown that Parts I and II are special cases of this problem. 
The basic approach will not be explained in detail here since it is given in Part I, but 
the problem will be solved independently of the earlier parts. 

1. Statement of problem. We deal with an incompressible fluid in a channel with 
a horizontal base, the fluid having been set in motion in such a way that the flow is 
two-dimensional and irrotational. The fluid motion is periodic with wave length \ in 


e) x 








> = > 














y 


Fic. 1. AB = }, OF =a, AD = h. 


the horizontal x-direction, and is assumed to be symmetrical about a vertical line 
through any crest or trough of the wave pattern. It is assumed that the wave profile 
travels from right to left (see Figure 1) with velocity c. The axes Ox, Oy move with the 
crest O, Ox being horizontal and Oy vertically downwards. Relative to this set of axes 
the fluid motion is steady. A velocity potential ¢ and a stream function y exist for this 
motion and are defined by 

d@ = udx + vay, 

(1.1) 
dy = —vdxr+udy, 


where (u, v) are the velocity components relative to O(2, y). 


*Received May 28, 1951. 
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From the symmetry assumption it follows that ¢ is constant along the verticals 
through O, B --- . We take @ = 0 along OF, ¢ = —@, along AD and @ = @¢ along BC. 
The lower boundary is a streamline which we take to be y = y, and the upper free 
surface is the streamline y = 0. Over the upper free surface the pressure p is everywhere 
constant. If g denotes the velocity of a fluid particle relative to O(z, y) and @ its inclina- 
: ‘", Suppose q is the relative velocity at O 


tion to the x axis, then u i = oO = ce” 


of the fluid particles, then Bernoulli’s equation applied to the free surface leads to 

Gg = go + 2gy, y = 0, (1.2) 
where g and y refer to values on the free surface. By differentiation with respect to @ 
Levi-Civita transforms this equation to the form 


00 qd re . 
—= — Hy sin 6, y = 0. 1.3) 


Oy c 


The method used here consists in replacing this condition by 


06 gl ieee 
. = = S54 sin 38, y = 0 (1.4) 
Oy c 
where / is a constant at our disposal. An alternative form of (1.4) is 
( ) 
J. dé gl 3¢\ ‘ 
Ch) — —e' 7 = @, y = 0 1.) 
| dw Cc ) 
where 9 denotes “‘imaginary part of”, w = @ + wand & = —7 + 76. 


now be stated as follows: £ is a regular 
in the @ direction, and is 


I 


The problem of the first approximation can 


function of win 0 < y < y, , is periodic and of period 2¢ 
such that 
a i{é} = 0 y=; 
b git} =0 , = 0 = +¢ 
s P QP Y (1.6) 
{dé gl | 
( j\2 + ae 7 =U Lb = (0 
dw Cc f ¥ 


In order to complete the statement of the problem it is necessary to postulate one 
further condition. In the infinitely deep channel problem and the solitary wave problem 
the actual fluid motion infinity is specified (¥ = © in the former, ¢ = +©@ in the 
latter) and in the present problem we must specify the actual flow at a convenient 


point of the fluid medium. This will be mentioned again later. 
2. The method of solution. We first write x = e ** and the problem in (1.6) may be 


restated as follows: x is a regular function of win 0 < y¥ < y, , is periodic and of period 
26, in the ¢ direction, and is such that 
(a) Jix} = 9, v=; 
a = _ 
(b) Six} = 0, ¢=0 od = +o 5x (2.1) 


i 5 a 3) 
(c) qe — (8g))/c i = 0, er 
x 
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Here, x is defined within and on the boundaries of an infinite strip and we consider 
next how we may continue x analytically throughout the whole of the w-plane. For 


this purpose we consider first of all the continuation of x across y = y, into the strip 
Yi < wy < 2y,. The line y = 2y, is the image of y = 0 in y = y, ; let P be any 
point of y 0 and let Q be the image point of P in y = y, . Suppose wp is the value 


of w at P and we that at Q. From Schwarz’s principle of symmetry, if we wish to con- 
tinue x analytically across y = y, , since d(x) = Oony = y, we must define the value 
of x at Y, namely x(weg), as follows 

x(We) = x(wp). (2.2) 
Here the bar denotes the conjugate, that is, if x(wp) = a + 78 then x(weg) = a — 7. 
The functions a and 6 vary with ¢ only as P moves along the line y = 0. The function 
x is defined in a similar way throughout the upper strip. In terms of @ and 8, condition 


2.1c) becomes 
da dg 3g!) 
B + —] = 0, = 0. 2.3 
' do + (22 t cc v \ ) 
It is easily verified that (2.3) is the condition also that the function 
[2 dx /dw (3ql)/c 
g(x) [2 dx/dw + (3g /e'| (2.4) 
x 
1as zero imaginary part on y = 2y, . Thus, if we write 
‘i [i dx/dw — (3ql)/e’ —— 
ju) = Ea (2.5) 
X 
then we have 9{f(x)} QOon Wy = Oand S{g(x)} = 0Oon py = 2y, , the only change 


being the reversal of the gravitational term. Thus we do not have perfect symmetry 
about ¥ yY, , but the result suggests, and we proceed tentatively from this point, 
that we build up the following image or periodic picture in the +y directions (see 
igure 2). 

In the breaking case a zero of x appears on y = 0 and we introduce, first of all, a zero 


of x at w |(a/t 3y,], where vy is a real parameter which is such that w/v > 3y, . 
We introduce a similar zero at the reflection of this pomt in y = y, , that is at w = 
ul (a/v y This we do because J(x) Qony = y, , in fact any distribution of 
singularities on @ 0 must be symmetrical about y = y, . We do not consider the 
image of the zero at w i|(r/v) — 3y,] in the free surface y = 0 from the analogy 
with the solution in Part I, but consider next the “image”’ of the zero w = i|[(r/v) — y] 
in the free surface y 0 where we must satisfy 9}/f(x)} = 0. The simple zero of x at 
w w/t ¥,| will, in general, be a simple pole of f(x) and we introduce, therefore, 
a simple pole at u —i|(r/v) — y,]. Similarly the “image” of the simple zero at w = 

i|(r v) — 3y,] in the line y = 2y, will be a simple pole at w = i|(x/v) + y,], and 
since this pole has to be reflected in y = y, , it follows that the poles at w = —7[(r/v) — ya] 
and w il(r/v) + y,] are double poles. Proceeding further, the image of the pole at 
w t{ (a/v ¥,] in y = 2, is a zero at w = i(x/v) + 3y,] while the pole at 
u m/v) + y,| will have as image a zero at w = —i[(r/v) + y,]. We now introduce 
it ¥ ¥, + 27/v and py = y, — 27/v lines along which 9(x) = 0 and at y = 27/» and 
y 2r/v + 2y, lines along which 9{f(x)} = O and g{g(x)} = 0 respectively. This 
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pattern is then repeated indefinitely in the +y directions along @ = 0. Similar patterns 


exist on all the lines ¢ = +2md¢, where m is any integer. There are no zeros nor poles 
of x at any other points in the w-plane. 

Since we have made g{f(x)} = 0 along y = 22/v and 9{g(x)} = 0 along y = 2, 
it follows that J(x) = 0 along y = y, + z/», this being the converse of Schwarz’s 
principle of symmetry. Hence g(x) = 0 on the perimeter of the rectangle FAHG and 





1{g(x)}-0 pray 
ye y+an 


\h° “ATTIILIINLILIILINIL 


k urei( (f +3) 
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x uw=-u(F -3¥,) 
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burs -i (B+) 
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Fic. 2. x indicates a zero of x, @ indicates a pole of x. 


this rectangle can, therefore, be mapped on the upper half of the y-plane and its perimeter 
into the real axis in the x-plane. The mapping function is Obtained from the Schwarz- 


Christoffel Theorem. Suppose the value of x at F is A; at G let x = A(1 — yw) and at 
H let x = A[l — (u/k)*], A, uw and k being real coihliies Then, since x is infinite at K 


the transformation required is the solution of the differential equation 
‘i / x \7? x \?/y Oe at 
ee ee ee . 
ax \ A \" A k - A : ead 


B being a complex constant in the general case. The appropriate solution of (2.6) is 


x = me = us| & (w — ty), «| (2.7) 
( od 
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where sn (a, k) is* the usual Jacobian Elliptic Function with parameter k, whose real 
period is 4K and whose imaginary period is 27K’. In obtaining (2.7) from (2.6) the 
constant B has been so chosen as to make x have period 2¢, in the @ direction; in this 


connection we use the relation sn (a + 2K) = — sna. Furthermore the arbitrary con- 
stant of integration arising from the integration of (2.6) has been chosen so that x is 
real when y = y, ; this enables us to satisfy (2.la). We can determine yu as follows. 


Since the poles of sn (a, k) occur at points congruent to iK’ and 2K + 7K’, on com- 
parison with Figure 2, we obtain 


i doK’ (2.8) 


v K 


We know also from Figure 2 that a zero of x occurs at w = i[(2/v) — y,] and hence 


i= a sn’| HK (: - 2v,), t| = 0 (2.9) 
ob Ww 


If we use (2.8) and the relation sn (a + 7K’) = 1/k sn a it follows from (2.9) that yu 


is given by 


Oo; Kk 
= k° an’ 2KY ’ t) (2.10) 


. ° . BE 37-310 
and thus the required solution is x = e “* = e”’ where 


— , 
x = {i — k’ an'( 2A ; k) a'| X (w — iy,), a} (2.11) 


3efore proceeding further we shall verify that x in (2.11) satisfies all the conditions 
enunciated in (2.1). The conditions of regularity in 0 < y < y, , of periodicity in the @ 
direction and condition (2.la) are certainly satisfied. To demonstrate that (2.1b) is 


satisfied we use Jacobi’s Imaginary Transformation, namely 


isn (u, K’) 
en (u, k’) 


sn (zu, k) 


, 


; 1 
, :) = ——_ 2. 
en (zu, k) mene” (2.12) 


se dn (u, k’) 
dn (iu, k) = —— 
mm en (u, k ) 


where k’ 1 — k*)'” is the complementary modulus. It then follows, when @ = 0, 
that x is real; similarly when @ = @ we have sn [(K/¢)(w — w,)] = sn 
(K + (iK/¢.)(W — ¥,;)] and since sn (a + K) = en a@/dn a the results (2.12) apply and 
make x real. Thus condition (2.1b) is satisfied on @ = 0 and @ = @¢, and elsewhere by 
periodicity. We consider finally condition (2.1¢c) or what is equivalent (1.4). 

The addition formula for the sn function leads to the result 


*When the parameter is not specified it will be understood that it is /. 
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62 
‘ K he sn (K¢/d,) en [iK(W — ,)/dbo] dn [7K(v — ¥,)/d0] 
sn — (w — iy) = ——— Se ; 
dbo 1 — k’ sn’ (K@/¢p) sn” [7K(~ — W)/do] 
(2.13) 
, en (K¢/do) dn (K¢/do) sn [7K(Y — ¥,)/d0] 
1 — k’ sn’ (K¢/q@) sn” [7K(~ — )/do] 
From (2.12) we see that en (7a) and dn (ia) are real while sn (7a) is imaginary, hence 
we obtain from (2.11 
a =e’ cos 380 
f 2iy, Kk \ sn (Ko/do) en” [i7K(Y — ¥1)/bo] dn* [7K(Y — yi)/do] 
= A4jl — k’sn } : 5 5 “ eer > 9 
( \ do / (1 — k' sn’ (Ko/q) sn” [7K (Y — ¥,)/do]] (2.14) 
en (Ko @Q dn Ko Do) sn” iK(yp = y,) doll 
1 — k’ sn’ (Ko/,) sn” [iK(W — y)/dol]? J’ 
8 = sin rf 
Gs K \ 
—e 2iy.K \ sate , —" T 
= 24k’ sn? | } sn (K@/¢,) en (K¢/d,) dn (Ko/do) sn [7K (CY — Y)/do] (2.15) 
\ @ 
en [tK(W — wv,)/dbo] dn [7K(~ — ¥,)/dbo] 
1 — k sn’ (K¢/¢,) sn [7K(W — y,)/bo]] 
Using (2.14) and (2.15) we can evaluate 06/dy. Since 7 tan 36 = B/a we have 
, 06 0D Oa : 
37e = = a — — B- (2.16) 
Oy Oy Oy 


and we substitute now for a and 8 from (2.14) and (2.15). To simplify the analysis we 


write 


? K¢ Ko K 
S =sn—, C= en D = dn a 
@ @ Do 
2.17) 
iKy iKy , iK y, 
8 sn : c, = cn : d, dn = 
rn Q Po 
We then obtain at y 0, 
3¢,c° (00/AaW)(1 — k’s?S 
2k*K A* sn® (2i¥,K/d)) SCD 
2iy,K)\,. 2 72 
= { | — }*S*g ke” s7 ( = ) Scd — ( D "4 f' cd, — 8d; _ k*sic; x 
" (2.18) 
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The detailed reduction of the right hand side of (2.18) can be omitted here and we 
quote merely the final simplified version of (2.18), namely 
_3doe"" (80/AY(1 — K*S*si)?_ _ (1 — 28; + h’si)(1 — 2k*st + ksh) | (2.19) 
2k*K A’ sn® [(2iy,K)/do]SCD — 1 — k’s} : 
If we now use the formulae for sn 2u, en 2u and dn 2u in terms of sn u, en u and dn u 
given in Whittaker and Watson’s Modern Analysis p. 498 Ex. 5 it follows that the 
right hand side of (2.19) is simply 


aie a “— 
(1 — k’s}) en (2K) dn (2m) 
Do Po 


so that 


06 2k°KA° 2 } { > 4 ] 1 : 
or 90|  _ 2RKA? (2K, — st) en (2K) 
oy y=0 o Qo do 


(2.20) 
di (24K) SCD 
d / (1 — ks S*)* 


From (2.15) we have 


——— — aK) SC Ds,c,d, 
e sin. = —2Ar me (—— —s 3503 ‘ 
a in 30 1k sn ( $y (1 — k’stS’y’ (2.21) 


y=0 


hence combining (2.20) and (2.21) we obtain 


a0 1 f2kK A en [(2ip, K)/$o] dn [(2iy,K) all - — 
MW lyo” 3 sn [(2i¥,K)/oo] e~sin 30. (2.22) 


If we now return to the fundamental boundary condition (1.11) we observe that it may 





be satisfied by (2.11) provided that 


3gl _— 2KA en [(2ipy, K)/do] dn [(27y, K) /g 
Cc) — oh sn [(27~, 4) /do] 


Bi (2.23) 


This is the formula for the velocity of propagation of the waves. Although the right 
hand side contains an 7 in the numerator it is real since en (7u) and dn (cu) are real 
while sn (7) is imaginary. 

This solution of the problem is completed by defining A. If we postulate that the 
actual fluid velocity at w = wy, is U, the relative velocity at this point is U + c hence 
from (2.11), since x = e** = (u — wv)*/c*, we obtain 


An (Uv + 1). (2.24) 


The solution of the problem can now be written in the form 


( w : , 2 r r r 
(#2) =(U+ on — k* sn° ( oe : t) sn° E (w — iy), |} (2.25) 


and the wave velocity formula in the form 


3ql 2iK en [(2ip,K)/do] dn [(27p, K) /do] 
__3gl__ __- 21K en [(21y,K)/4o] dn bo] 2.26 
(U +0) po sn [(27y,K)/$o] satiate 








64 T. V. DAVIES [Vol. X, No. 1 


This wave velocity formula can be expressed in fairly simple terms by using the Jacobi 
and Landen Transformations. In the first place, applying the Jacobi transformation to 
the right hand side of (2.26) we obtain 


3ql 2K dn [((2y,K) /do, k’] ae (2.96’) 
(U + o> sn [(2¥,K)/do, k'] en [(2¥,K)/do, k’] ~~ 
Landen’s transformation states that 
f i 1— k’ I he sn (u, k) en (u, k) 
Sli + & )U, —sT( = = igh 3S eee = ——s. 
l+k dn (uw, k) ‘ 
hence 
2 >] = . PY Aa 
= : , 20) a k)K Po _ oe = a (2.27) 
U+e sn {2(1 + k)Ky/b , (1 — &)/1 + &)} 


If we multiply both sides of this by y¥, the right hand side is then in non-dimensional 
form. 


The breaking case is at once evident from Figure 2, for this corresponds to the case 


when the zero at w i[(r/v) — 3y,] lies on y = O, that is the breaking case occurs 
when 
T 
= SY, (2.28) 
y 


Using 2.8) this may be written in the form 


Ky, eee 

a : % 99 

= 5K’. 2.29) 
Po o 

This result is easily verified also by writing x = 0 in (2.11) and solving for Ay,/¢o 

when w = 0. This quantity 7/y is one of the important parameters of the problem for 


we may note also that the small amplitude problem corresponds to large values of 2/v, 
when the system of singularities moves away from the region 0 < y < y, 

There are two points which must be mentioned before concluding this section, the 
first is that the problem of the uniqueness of the solution (2.11) has not been solved. 
Secondly, the problem of higher approximations has not been attempted, but it seems 
fairly certain that the higher approximations may be derived as in Part II by using 
the more exact boundary condition 


= g {1 a 5 30 + =e" sin’ 30 + ao e-™ sin® 30 + | (2.30 
ay — VET ere ; = 


and to assume a solution of the form 
= f K ' 6| AK ‘ 
¢ . = Ay] T MM, SN i ry) + M3 Sn Tes ry) + hn i 2.31) 
; Po 


where », and pu; are constants. The results of Parts I and II suggest that —y, is the same 


as the uw of (2.10) and 
1 « /2K 
vs = ~-,5n ( vs ; k ; 
54 Po 
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also if u, and higher terms are neglected, the results are likely to be in error by about 
13% at the breaking point but will be much less for a non-breaking wave (here we 
assume that / 1/3). 

An alternative method of deriving the solution (2.11) is to make use of the isolated 
wave singularity distribution, given in Part IT, and to use theta functions for the more 
extended distribution of the present periodic wave problem. This has been done, but 

no shorter than the method expounded here. 

3. Special cases of the solution (2.11). We demonstrate first that the problems of 
the infinitely deep channel and the solitary wave are merely particular cases of the solu- 
tion (2.11) to the order of approximation in this paper. For this purpose we use the 
hetween the Jacobian Elliptic Functions and trigonometric functions, 


connections he 


namely 
(2Kz , 2r fq smxv , q *sin 32. 
T / Kk \ 1 q l-—¢@ 
9K >» = " » . .... 
en 2Ka 1) Zi q COS 2 L q COS 8a L 
\ Kki l+q — a ow i 
/ 
, (2K2 t) 2 27 Jqg cos 2x 4 ¥ 0S 4x +} 
\ nr? 2K Kil+q 1+¢@ 
1 — k’ 2y + 2q° + 2q° + 
] , £ Ss € l | 
| / 1 + 29° + 2q'" 4 
\ end » Zz ends to zero and we have in the limit 
; i; 
lim ( i73j = | (3.1) 
o 4 
furthermore as 0, K — x/2 and Kh’ — >=, hence 
/ , » 
(2h ; 2K2 ; 2A x ’ 
sn | }=sinz, lim en ( ,k}) = cosx, limdn ( ,»k}=1. (8.2) 
Let us consider no hat happens to the solution of the previous section as the parameter 
tends to zero. When ¢, is constant we see from (2.8) that z,v tends to infinity, that 
is, the system of singularities moves to an infinite distance from the stripO < ~y <p. 
In the problem of Part I, namely y, , there is a line of zeros parallel to y = O lying 
in the half plane ¥ 0 and to reproduce the results of Part I we, therefore, make k — 0 
nd ¥, © in such a way that one line of zeros remains in the finite part of the w-plane. 


As / > 0 we thus obtain 


K wm ° \ I ( rw/2o Tyi/2¢ —trw ry °) 
sf > — ¢ ( i 
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hence as k — 0 and y 


Thus from (2.11), ask — 0, y,; - 


~ Ad = (4 ke iste’ Joi 3 } (3.3) 


The two limiting processes are independent of one another; if we now postulate that 


: er 
lim lim — k’e’”’ = constant (3.4) 
> 0 ) 


and make U’ = 0 in (2.24), we obtain the solution of Part I, namely 


e"* = |] — Be**”’. (3.5) 
The velocity of propagation formula (2.26) may be discussed in a similar way. With 
U = Othe left hand side is 3gl/c° and the right hand becomes 
ir cos [(r2y,)/do] 


Po 
do sin | (rz) do | 


as k -0, so that as y >co we obtain 


3gl ul — 
——- . — (3.0) 
4 @D 


which is similar to Part I. 
Consider next what happens when we make k tend to its upper limit, namely unity. 
As k — 1, k’ — 0, and thus we can deduce the relevant limits from Jacobi’s Imaginary 


Transformation. We obtain from (2.8) and (2.12): 





[K ) (vK’ 
sn 4 w— tW,), ke = sn‘ (w — ty), k 
\d ) t T 
.sn }(vK’/ir)(w — ty,), k’} 
a ae anal ' 
en {(vK’/ir)(w — ty,), k’} 
and as k’ — 0, since K’ — (1/2)z, we obtain 
K ; | l i 
lim sn {— (w — ty,), ke = tanh 5 ¥(w — ty). (3.7) 
; \do ) 2 
Similarly, 
2iy,K ; 
lim sn ( “ , kk} = ttanvy,. (3.8) 
—1 0 


It will be noted that by first substituting for the ratio K/¢, in the above expressions 
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we avoid the difficulty of indeterminacy of this ratio. The solution (2.11) now becomes 


x = Afl + tan’ vy, tanh’ 4o(w — 7y,)} 
(3.9) 


= A sec’ vy,{1 — sin’ vy, sech’ 4x(w — iy,)} 


which is the isolated wave solution of Part II, A being equal to cos’ wJ, . Similarly the 


wave velocity formula (2.26’) becomes 








Bgl (2 cos vv) [(K’y) ‘| 7 
c* sin [(2¥,vK’)/x] cos [(2ywK’) r|’ 
hence 
3g! = vp cot ry, 
as In Part II. 


We can summarize the results of varying parameters as follows: 

(a) k — 0, ¥, ©; finite amplitude waves in a channel of infinite depth; 
b) | > 1: isolated wave: 
(c) k — 0; ¥, , » finite; small amplitude waves in a channel of finite depth; 
d) Ay, @ 1K’; breaking condition at the crest. 
e) 0 Ky,/K'o < 1/3. 
The result (e) follows from expressing the condition that z/v must lie between 3y, and 
infinity. The result (c) follows from the early part of this section for in this case formula 
(2.11) becomes 

y 

t i » - 9 ° 
et = Ag + k’ sinh’? — sin’ > (w — iv) 
\ bo 2, 


and, k being small, this approximates to the usual solution for small amplitude waves. 
Furthermore the wave velocity formula in this case becomes 


q T T 
= coth va 
c do 0 
and with the usual assumptions ¥, = ch, @ = $cd this becomes 


2r 


c= gh tanh (24) 
y= ‘ we 

















—NOTES— 


ON SOME STATISTICAL PROPERTIES OF HYDRODYNAMICAL AND 
MAGNETO-HYDRODYNAMICAL FIELDS* 


By T. D. LEE (University of California, Berkeley) 


1. Introduction. It is now recognized through the works of Taylor, Kolmogoroff, 
Heisenberg, Batchelor and many others that turbulence, rather than laminar flow, is 
statistically the most natural phenomenon exhibited by a non-conducting fluid. From the 
knowledge of magneto-hydrodynamical waves,’ it is also believed that in a fluid of high 
electrical conductivity, magneto-turbulence may be the most natural phenomenon. The 
purpose of this paper is to study some of the statistical properties of a fluid and their 
relations to the theory of turbulence. The discussion of magneto-turbulence leads to the 
conclusion that in the usual Kolmogoroff region, where the behavior of the fluid is inde- 
pendent of detailed mechanisms of energy supply, the magnetic energy will be in equi- 
partition with the kinetic energy and that the energy spectrum of either field will be 
proportional to k *’” where k is the wave number. 

The ergotic motion of an idealized fluid with no dissipation of heat is studied in 
section 2 and the existence of a Lionville Theorem is proved. The application of a dimen- 
sional consideration and some of these statistical properties to the steady state of magneto- 
turbulence is discussed in section 3. 

2. Ergodic motion of an idealized fluid. Let us first consider an incompressible fluid 
with zero viscosity and zero electrical conductivity. If a certain amount of energy is given 
to this fluid, we wish to know statistically how this energy will be subdivided into various 
modes of possible motion. Let v(r), p(r), p be the velocity, pressure and density of the 
fluid. The Euler’s equation of motion givest 


v = -(v, Vv — ~2 (1) 


and 
V-v=0. (2) 


It is convenient to analyze the motion of the fluid in terms of the Fourier components 


of its velocity. We write 


p(t) = ple" (3) 
v(t) = DL v(de™* = YO [a(ke) + a9(le)Je** (4) 
with p(k) = p(—k)*, v(k) = v(—k)*, a(k) = a(—k) 
and g(k) = —g(—k). 


*Received April 2, 1951. 
1H. Alfven, Arkiv. f. Mat. Astr. o. Fys. 29B, 2 (1943). 
tA time derivative is indicated by a dot following the letter. 
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Equations (1) and (2) can then be written as 


p(k) = — 7 2 [v(k + &’)-k][v(—k’)- ki], (5) 
v k’ 
v'(k) = —i D [v(ke + ke) kl(—k’) + 28D [v(e +k’) K]fv(— ke’), (6) 
k’ = = 
and v(k)-k = 0. (7) 


For mathematical convenience we shall treat the three components of a(k) and 6(k) as 
independent but regard (7) as a constraint applied to the initial condition of the fluid. 
For from (6), if (7) is true at a particular moment it is always true at any other time. 
Let us now consider a phase space with a,(k), a,(k), a,(k), 8,(k), B,(k), 8.(k). ---, 
a,(k’), a,(k’), --- . as its coordinate axes.” In this space each point, compatible with the 
initial condition (7), represents a dynamical state of the fluid. The trajectory of this point 
governed by (6) describes the subsequent motion of the fluid. Differentiating (6), we have 


Bai(Ik) , 98;(kk) _ 


0; 6 = 2, ¥, 2. 8) 
Oa ;(k) 08 ;(k) 8) 


This means that the density of a group of representative points in this phase space remain 
constant along their trajectories. Just as in an ordinary mechanical system, the existence 
of this Lionville Theorem suggests that the a priori probability of finding this system in 
any region of this phase space is proportional to the volume of this region.* From (6) the 


law of conservation of energy can be written as 


b 3 [a°(k) + 8°(k)] = constant (9) 
k 


Thus, we find that statistically every mode of the Fourier components of velocity must be 
in energy equipartition. Let F(k) be the spectrum of kinetic energy per unit volume. In 
the equilibrium case, we find 
7 2 
F(k) « k’. (10) 
Next we discuss the corresponding case of an incompressible fluid with zero viscosity 


and infinite conductivity. Let j(r), E(r) and H(r) be the electric current, electric field and 
magnetic field respectively. The hydrodynamical equation becomes 


, Vv ixH 
v= —(7, Vw — 2+ tAS. (11) 
p pe 
Due to the infinite conductivity, we have 
E=-"xH. (12) 


2To avoid the difficulties of an infinite number of degrees of freedom, we may introduce in (3) and (4) 
an upper limit on & corresponding to, say, the inverse of atomic dimensions. 

sAlthough the ergodic hypothesis may be true for a three dimensional fluid it is not true for a two- 
dimensional case. The existence of a conservation law on vorticity forbids the fulfillment of ergodic 
hypothesis for a two-dimensional fluid. See T. D. Lee, J. App. Phy. (in press). 
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Furthermore, j, E and H are connected by the usual Maxwell Equations. On eliminating 
j and E from (11), (12) and the Maxwell Equations, the equation of motion of this fluid 


can be written as 


' Vv H’\ 1 
v = —(v, Vw — — (p+ 2) + 7 @, *)E, (13) 
p Sr 4rp 
H' = —(v, V)H + (H, V)v (14) 
and V-y= V-H = 0. (15) 


The term 
2 
-v(@)4 1a, wH 
Sr Ar 


in (13) is the force acting on the fluid due to Maxwell stress. In an infinitely conducting 
fluid the magnetic lines are always frozen in the matter. Equation (14) means that in a 
system moving with the matter the increment of H is given by the extension of the mag- 
netic lines of forces due to relative motions of matter in different places. 

As before, we make a Fourier expansion of the field variables. In addition to (3) and 


(4) we have 


H = >> H(kje™* = D> [g(k) + th(k)]e™" (16) 
k k 


with g(k) = g(—k) and h(k) = —h(—k). 
In terms of these Fourier components (13), (14), (15) can be written as 


vik) =i>d {—Iv(k + k’)-k]w(—k’) + = (H(k + K’)-ea(— «')} 


(17) 
1 ( 
+45 5 lw(e + i’) e}fv(— ke’) om [H(k + ke’) = eEK(—’)- eh 


H'(k) = i D0 {—[v(k + k’)- kJH(—k’) + [H(k + k’)-k]v(—k’)} (18) 


v(k)-k = H(k)-k = 0 (19) 


and 
From (17)-(19) it can be verified readily that 


dai(k) , 98;() 
da,(k)  98,(k) 


= 0, 


and 
dgi(k) , dh; (ke) 
ag) ah, (kk) 


= Q; t= 2, y, 2. (20) 


The law of conservation of energy now becomes 


by . [a’(k) + 6°(k)] + = lg’(k) + h'(19)} = constant. (21) 


k 
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Let us now consider the phase space with a,(k), 8,(k), g.(k), h.(k), a,(k), 8,(k), g,(k), 
h,(k), «+ , a,(k’), --- . ete. as its coordinate axes. Every point in this space, compatible 
with (19), represents a dynamical state of this infinitely conducting fluid. Equation (20) 
suggests again that the a priori probability of finding this system in any region in this 
phase space is proportional to the volume of that region. Hence we find that in the 
equilibrium distribution every mode of the Fourier components of magnetic field and 
velocity field must be in energy equipartition. Let 1/(k) be the corresponding energy 


spectrum of magnetic field per unit volume, then we have 
Mik F(k) « k’. 22) 


3. Turbulence and magneto-turbulence. In the case of a real fluid, due to the energy 
dissipation by either vist osity or electric resistance, a steady state is reached only if there 
is a constant source of energy supply. This situation can best be represented by a con- 
stant energy flow from small wave number regions to large wave number regions and then 


in turn a transformation into heat. The existence of this energy transfer makes it neces- 


sary that the energy spectrum must be quite different from the equilibrium case of an 





idealized fluid. However, since this energy transfer is only between different | k | values, 


it is expected that all modes of Fourier components of the same | k | values may be still 


in energy equipartition. We shall discuss only the statistical properties of magneto- 
turbulence in some detail as the corresponding properties of turbulence are well known 


and also can be obtained readily from the following discussion by putting H = 0 


A. Steady case of nviseid infinitely conducting fluid with a constant energy flow. Let us 
consider again an incompressible fluid with zero viscosity and infinite conductivity. If an 
amount of energy is supplied at a certain wave number, then due to the existence of 


non-linear terms in the equations (17) and (18), this energy is scattered between the 
magnetic field and the velocity field, and also is transferred within the different wave 
numbers of both fields. Due to the availability of more volume in phase space for high & 


values, statistically this ene rgy is always transported toward higher and highe wave 


numbers. Let us now impose a boundary condition that a constant rate of energy ¢ per 
unit volume is supplied at wave numbers smaller than a certain k, and the same rate of 
energy is taken out at I, much larger than k, . Since the transfer of energy from low k to 


high & is of a statistical nature, it is expected that there exists a significant region, defined 


by limits ko and ke (ko < ke) between k, and k, such that in this region the behavior of the 


field depends only on ¢ but not on the boundary and the detailed mechanism of energy 
supply. This region will be called as ‘“‘Kolmogoroff region” since Kolmogoroff* was the 
first one who pointed out the existence of such a region in ordinary turbulence. Let 
E,(k) and Ey(k) be the rate of energy tr: 


g (I sfer per unit volume at k, within the velocity 
field and within the magnetic field respectively. In the Kolmogoroff region all functions 
| ; } 





can depend only on ¢, p and k&. From dimensional considerations the only forms that 
F(k), M(k), E,(k) and Ey(k) can take are 


F(k) ee p’ Pa 37.—5/ A (23) 
M(k) = Cyp'2h-”, (24) 


E(k) = «,€, (25) 
‘A N. Kolmogoroff, A.N.C.R. Acad. Sci. U.R.S.S. 30, 301 (1941); 32, 16 (1941). 
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E(k) = kué, (26) 


with Ky =f Ky = Es (27) 


C,, Cu, k. , Ky are numerical constants. Equations (25) and (26) mean that the rate 
of energy transfer /,(k) and Ey(k) must be constant at every k. Consequently, there 
can not be any energy transfer between the magnetic energy and kinetic energy in this 
Kolmogoroff region. From the result of previous section we have that for the equilibrium 
distribution every mode of the Fourier components of this idealized fluid must be in 
energy equipartition. A deviation from this distribution necessarily requires an energy 
flow. Hence, the presence of energy flow from low | k | to high | k | makes the spectra F'(k) 
and \/(k) different from k° law. Similarly, since there is no energy transfer between the 
magnetic field and the velocity field we can conclude that for any fixed value of | k | the 
magne energy must be in equipartition with the kinetic energy. We may write 

C, = Cz. (28) 


Similar arguments also lead to the conclusion that both magnetic field and velocity 


field are isotropic in the Kolmogoroff region. If the lower limit of Kolmogoroff region ko, 
is comparable to k, and the upper limit ke, is much greater than ko, one may normalize 


F(k) in terms of the average velocity, v) , of the fluid. We have 
gu | ik) dk. 
Hence, 
M(k) = F(k) = ; pk *vsk-*”* (29) 


in the IKXolmogoroff region. 

B. Effects of viscosity and conductivity. The existence of viscosity and finite conductivity 
gives a particular mechanism of energy dissipation in the above discussion. Let ¢ and o be 
the viscosity and electric conductivity respectively. We have the expression for energy 


dissipation per unit volume, [€]aissipation - 








. s2 ss ,2 | Vv H 12 
[€]aissina =y|V xvPp+hep vVxvli Ch tele (30) 
o (4ir)"o 
Equation (30) can be written in terms of the energy spectrum as 
"(k Quk* 
EQ =—-~ Fk) (31) 
Gt Soisciier p 
and 
aM | ck? 
—- = — — M(k). 32 
| dt finite conductivity 2ro ( ) ( ) 


The dependence on k’ in (31) and (32) means that the effects of viscosity and finite con- 
ductivity are important only at high k values. The upper limit ke of the Kolmogoroff 
region may be estimated as the wave number at which the energy transfer between wave 
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numbers is comparable to the heat dissipation. From (17) and (18) the efficiency of energy 
transfer between different wave numbers can be estimated as 


- : ] dF (k) ae ey ; l : dM (k) ai TRL. X 
Ky Ee dt 7 _ Ky ais dt i a vi “e iP (33) 


Hence kc is given approximately by 


(4 +) ~ VFI. 4) 


Using (29), (34) can be written as 


k 


; = R* (35) 


where F is the modified Reynolds number defined by 


R= an / (u + “2 \i, . (36) 
\ 4ro 


Equation (35) means that the range of the Kolmogoroff region increases with increasing 
Reynolds number. The behavior of the energy spectrum between ko, the lower limit of 
the Kolmogoroff region, and the boundary k, must depend on some extra-dimensional 
quantities related to the boundary or the detailed mechanism of energy supply. Hence no 
universal energy spectrum can be given in that region. The importance of this transition 
region will depend on the efficiency of energy coupling between the magnetic field and the 
velocity field and also on the detailed manner as how the energy is supplied. In ordinary 
turbulence it is assumed that such a transition is unimportant and this assumption is 
experimentally verified.” It may be natural to make the same assumption here. Due to 
the lack of experimental evidences, such assumption is always subject to criticisms and 
future verifications. But if this is the true case, then at high Reynolds number the total 
mechanical energy must be comparable to the total magnetic energy in a conducting fluid. 

The author is indebted to Professors E. Fermi and 8. Chandresekhar for discussions. 


5See e.g. Heisenberg, Zeit. f. Phy. 124, 628 (1947). 


ON THE INVERSION OF THE VOLTERRA INTEGRAL EQUATION* 


By B. GROSS (The National Institute of Technology, Rio de Janeiro) 


When the given kernel of the Volterra integral equation can be represented as a 
Laplace transform, the same representation is obtained for the resolving kernel of the 
equation. For this case the solution is given in explicit form. 

The Volterra integral equation 


az 


fz) = gz) + | Ke-wow dy @>O) (1) 


0 





*Received May 8, 1951. 
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has proved to be of interest for the solution of many problems in physics. Its solution 
can be reduced to that of the simpler equation’ 


: 
K(x) = mx) + [ ke — dm(O at. (2) 
Jo 
Here, k(x) is the given kernel and m(zx) the resolving kernel of Eq. (1). Volterra’ has 
given a solution of this equation in form of a series development. Doetsch* has solved 
it with the aid of the Laplace transform and Titchmarsh*® with the aid of Fourier trans- 
forms. All these solutions however allow for very general conditions to be imposed 
upon the given kernel. 

In physies integral equations of this type mainly occur in the theory of dielectric* 
and elastic” relaxation processes. In these applications k(x) belongs to a very particular 
class of functions; it can be shown to satisfy conditions which will allow it to be repre- 
sented as a Laplace transform 


a oO 


k(x) = | r(s)e"*" ds. (3) 

The purpose of the present paper is the discussion of the solution of Eq. (2) under the 
condition that k(x) is given by the expression (3). 

Equation (2) is of the convolution type. Application of the Laplace integral method 


of inversion’ gives therefore 


M@) = 7 or (4) 
where 
K(p) = [ e"*k(z) da, (5a) 
M(p) = [ e”*m(x) dx. (5b) 
It follows from (3) and (5) that 
K(p) = [ e”* dx i" r(s)e"*" ds. (6) 


Suppose that r(s) is a real function of integrable square over (0, ~). Then the order of 
integration in Eq. (6) can be inverted and 


K@ = | =. ds. (7) 


1V. Volterra and S. Pérés. Théorie des fonctionelles, Paris, Gauthier Villars, 1938. 

2G. Doetsch, Math. Annalen 89, 192 (1923). 

8E. C. Titchmarsh, Introduction to the theory of Fourier Integrals, Oxford, Clarendon Press, 1937, 
p. 311. 

‘B. Gross, Zeits. f. Physik 107, 217 (1937). Phys. Rev. 57, 57 (1940). 

5B. Gross, J. Applied Phys. 18, 212 (1947); 19, 257 (1948). 

6R. V. Churchill, Modern operational mathematics in engineering. McGraw-Hill, New York, 1944, 


pg. 42. 
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K(p) is therefore a Stieltjes transform. As such it is regular in the p plane cut along 
the negative real axis and the functions K(pe™*") are also of integrable square. 
Applying the complex inversion theorem to Eq. (4) one obtains 


P)__ ops dp, (c > 0). (8) 


Suppose also that 

1+ K(p)|>6>0. 
Then the line of integration can be deformed into a loop round the negative real axis, 
so that ultimately 


ee E= )_ _ _Kige | dp. (9) 
Jo 1 + K(pe *") 1 + K(pe"**) 


att 


Equation (9) can be written in the form 


a oO 


m(z) = | r(pje”* dp. (10) 


This shows that under the given conditions the resolving kernel too is represented by 
a Laplace transform, and 
I K(pe'*”) K(pe™*") (11) 
= = area roe \ 
2rt L1 + K(pe'*) 1+ K(pe**”) 
Here, K(pe™'*) is a shorthand way for writing lim;.. K(—p + 76). Then, according 


to Eq. (7) 


rrr ‘ ' r(s) 
K (pe )= lim | a ws i8 ds, (12) 
and 
a 4 ‘(g) a *(s) \ . 
lim | e = #G 16 ds = | ar D ds - inr(p), (13) 


where the integral denotes the principal value.‘ Substitution of (13) into (11) and 
decomposition into real and imaginary parts gives 


o “( 
rp) = ; rp) = — (14) 
[ ") oo. as| + 2°*(p) 


Equations (3), (10), and (14) represent the solution of the given equation. 
The author is very indebted to Professor E. C. Titchmarsh whose help and advice 
has enabled him to write the present paper. 


7The evaluation of this limit represents a particular case of the calculus of the boundary values of a 
J fp 

function along a given path. The expression (13) is immediately obtained from the general formulae 

given e.g. by A. Hurwitz and R. Courant, Funktionentheorie, Julius Springer, Berlin 1929, p. 333. 
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ON CRACKS AND DISLOCATIONS IN SHAFTS UNDER TORSION* 
By ALEXANDER WEINSTEIN (University of Maryland) 


1. Introduction. Several decades ago W. Arndt [1] showed that the theory of shafts 
of revolution under torsion can be reduced to the investigation of axially symmetric 
motions of a fictitious incompressible fluid in a space of five dimensions. Let (x, y) be 
the meridian plane of the shaft, x being the axis of symmetry. All relevant quantities 
are determined by the potential @(z, y) and the stream function (zx, y) which satisfy 
the generalized Stokes-Beltrami equations 


yb, = V,, yb, = —W,. (1.1) 


The displacement of a particle is perpendicular to the meridian plane and is given by 
the formula v = y®(z, y). The stresses in the meridian plane are rs, = y®, and 7», = y®,, 
where @ denotes the polar angle in the plane perpendicular to the z-axis.** Any stream- 
line V const. is free from stress and can be taken as the profile of a free surface 
bounding the shaft. 

By analogy with ordinary hydrodynamics, Arndt developed the method of sources 
and sinks in a space of five dimensions using however source distribution only on the 
axis of symmetry. Recently the method of sources and sinks has been extended by the 
present author to distributions outside of the axis of symmetry [2], [3]. This work was 
extended by A. Van Tuyl [4] and -by M. A. Sadowsky and E. Sternberg [5]. The present 
paper gives new application of this theory to problems of cracks and dislocations in a 
shaft under torsion. A new and simple correspondence principle will be used connecting 
the method of sources and sinks with electrostatic problems in Generalized Axially 
Symmetric Potential Theory [6], [7] which will be denoted by the abbreviation GASPT. 

2. The correspondence principle. GASPT deals with equations of the form 


y’¢(z,y) = ¥,(2,y), y’e,(z,y) = —¥.(2, y) (2.1) 


for all non-negative values of the parameter p. For p = 1, 2, 3, --- ¢ represents an axially 
symmetric potential in the meridian plane of a space of p + 2 dimensions. The case 
p = 0 corresponds to harmonic functions in the plane. In all cases y is the corresponding 
stream function. For clarity’s sake we shall sometimes write ¢{p} and y{p} in place 
of ¢ and y. We have from (2.1) the following equations 


J » de) a ( » 9e) a ( le 14 a ( -» a) 
4+ — —)]=0, — = — 7—-} = 0. 2.2 
Ox (y Ox Oy y Oy Ox y Ox + Oy y Oy (2.2) 


We see that ¢ and y satisfy the same equation only in the case p = 0, which corresponds 
to the well known fact that in plane hydrodynamics a stream function can be used as 
a velocity potential and vice versa. For p = 0, we use a fundamental identity [6, p. 351, 


(48)], which can be written as follows 
vip} = Cy’*'elp + 2}, (2.3) 


where C is an arbitrary constant. This identity permits us to obtain from a stream 
function ¥{p} a certain potential ¢{p + 2} and vice versa. 

*Received June 7, 1951. This work has been done under the sponsorship of the Office of Naval 
Research. 

**We have put Lamé’s constant » equal to unity. 
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3. The method of sources and sinks. In the case of shafts, as well as in any other 
cases, solutions of certain problems are obtained by putting, for a constant U, 


{ip} = Ux — ofp}, Vip} = U(p+1)'y’"" —v{p}. (3.1) 

In the usual cases ¢ and y are due to sources distributed in a finite domain. For the 

theory of shafts, p is to be taken equal to 3. Our main problem is to determine the lines 

V{p} = constant, the line V{p} = O being the line enclosing the singularities. We 

assume in the following that this line is closed. By (2.3), we may put y{p} = 

: aah ne rp ; : 

U(p + 1) ‘yelp + 2}. Using the notation ¢{p + 2} = ¢*(z, y), we have (3.1) 

Vip} = U(p+ 1) 'y’"'01 — ¢*(z, y)). (3.2) 

This shows that fo the closed streamline V{p} = O correspond the level line ¢* = 

¢{p + 2} = 1. As ¥{p} tends to U(p + 1)"'y’”’ at infinity, we see from (3.2) that 


¢g{p + 2} must tend to zero for 2° + y* —~&. In this way the problem for ¥{p} is now 
reduced by our correspondence principle to the electrostatic problem for g{p + 2}. 

4. The crack problem. A flat interior circular crack perpendicular to the axis of 
symmetry of the shaft is defined as a circular disc separating the material. It is assumed 
that no stress is acting across the surface of the disc. The definition of the crack adopted - 
here is the same as that used in another problem by M. A. Sadowsky and E. Sternberg 
[8]. Let us note that A. E. Green [9] uses the word crack in a different sense, namely in 
order to denote the displacement on the surface of a body due to localized pressure. 

Let x = 0, y = 0 be the center of our disc of radius b. By Arndt’s theory, the boundary 
condition is given by ¥{3} = 0 on both sides of the segment « = 0,0 S y S bin the 
upper half of the meridian plane, ¥ being an even function of z. The problem is, by 
(3.2), reduced to the determination of the electrostatic potential ¢{5} of a disc of radius b 
in a space of seven dimensions. This problem will be now solved by using ellipsoidal 


coordinates. Let 


9 


x ay 
sl + “= Y 7 = l (4.1) 
rn A+ b 
be a family of confocal ellipses which generate through a rotation about the z-axis a 
family of oblate spheroids \ = constant in a space of seven dimensions. The value \ = 0 
corresponds to the disc. Let us put 

« _ 80° [ du (4.2) 

: 39 Jy (u+ 0B°)*(u)'” ; 
Then ¢* satisfies (2.2) for p = 5, vanishes at infinity and takes the constant value one 
on the disc. The integral in (4.2) can be expressed in terms of elementary functions, 
the result being the formula 
‘asks 2b (: 5b’)(a)'”? 
g* = _—— Q a — Me a (4.3) 
T b 3r (A + Db) 
By Sec. 3 we have 
y{3} = 4°°Uy'*¢* (4.4) 

In order to obtain ¢{3} we introduce in the z, y-plane elliptic coordinates by the formulae 


x = bsinh é cos 7, y = b cosh é sin 7. (4.5) 
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By comparison with (4.1), see that 

2: 192 2 2 2 

b° sinh” — = 4, b° cosh” —§ = A+ 6 (4.6) 
As (4.5) defines a conformal transformation, i.e. a local rotation and uniform dilation 
of the coordinate system, we have, in place of (2.1) 


yee =y¥,,; Y' Pr = — ’ (4.6) 
where by (4.5), y = 6 cosh é sin 7. 
By (4.3) and (4.4) we have 
¥{3} = Ub‘(6r)~'{3 cosh‘ £ arccot (sinh ¢) — 3 sinh® — — 5sinh €} sin‘ 7. (4.7) 


To find ¢, we use (4.6) for p = 3, which yields ¢, . By an integration we obtain finally 


o{3} = Ue + 2sinh é arccot (sinh ~) — 2} COS 7. (4.8) 


& 2 
cosh’ £ 


OO 


Obviously ¢ vanishes for x + y° —@. 
The formulae (4.8) and (4.7) together with (3.1), taken for p = 3 give the solution 
of the crack problem in elliptic coordinates. The same solution can be also expressed 


in terms of x and y. In fact, the quantity 


s(x, y) = [ae tly—bPf?+iet+tiyt+o Py” (4.9) 
is a function f(A) of \ alone. Putting y = 0 in (4.1) we find \ = 2°. Therefore on the 
ellipse \ = constant we have s = s(x, 0) = 2(x° + b)' 7 = 2(\ + b’)'”. This equation 
yields 


h = 43s°*(2, y) — B’. (4.10) 
This value of \ can be substituted in (4.6), (4.7), and (4.8) and gives us explicit formulas 
for g(a, y) and W(x, y) in terms of x and y. 

The same functions ¢ and y can also be expressed in terms of integrals involving 
products of Bessel functions. However, such integrals are less suitable for the discussion. 
In fact, these new expressions should be rather considered as new formulae yielding 
certain integrals in terms of elementary functions. The corresponding results (also for 
an arbitrary value of the parameter p) will be published elsewhere. 

5. Discussion of the results in the crack problem. For the discussion, we have to 
use the formulae for displacements and stress given in the Introduction, and the formulas 
for &, Y, ¢ and y given in paragraphs 3 and 4. It will suffice to note the following results. 


The transverse displacement v(x, y) is given for x = 0 by the formulae 


jad Uy(b? — y’)'”, for 0Osy<sb, 
e T 


0, y) = (5.1) 

l 0, for y> ob. 
The signs + and — correspond here and in the following to the positive and to the 
negative side of the y-axis respectively. The stress 79, is zero on the disc, i.e. for x = 0, 


0 < y S b but becomes infinite like (y — b)~'”” for y > b. The stress 7», tends to infinity 
in the same manner for « = 0, y 7 b. For large values of x” + y’, ¢* tends to zero like 


(2° + y°)°”. Therefore the formula ¥V = 4°'Uy‘(1 — ¢*) shows that already for 
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moderate values of y, a streamline Y = constant is essentially parallel to the axis of 
symmetry, so that we have a nearly cylindrical beam with a circular crack. 

It is clear that the method presented here can be immediately applied to cavities 
of the form of oblate spheroids. 

6. Dislocations. According to the classical theory of Volterra a surface of dislocation 
is a surface of discontinuity for the displacements. The stresses remain continuous 
across this surface and the displacements on both sides of a dislocation differ only by a 
rigid body motion. In the case of a shaft under torsion we take as the surface of the 
dislocation the same circular dise as in the crack problem. The only admissable disloca- 
tion is, in this case, a rigid rotation of the negative side of the disc relative to its positive 
side. Such a dislocation is readily obtained by taking in (3.1) for ¢ = ¢{3} the potential 
of a uniform magnetic shell in a five dimensional space. By a specialization of a general 
formula [7, p. 76, (6.1)] we have 


og = Mb’'y" | e*Ji(yt)J.(bé) dt, c=0 (6.1) 
and 
y= —Mby? | c*'Ji(yt)J(bt) dt, 2x20 (6.2) 
For x < 0, we have ¢(—2, y) = —¢(2, y), W(x, y) = ¥(—2, y). In (6.1) and (6.2) M is 
a positive constant and the J, denote Bessel functions. In view of the results given in 
[7] we have the following values for the displacement v in the plane of the dise: 
yey, for y<d 
0, y) = < (6.3) 
( 0, for y > b 


Moreover we have 
lim v(a, y) = Mb(r — 8)x™* 

for zx — 0,7 — band (b — y)x”' — cot 8, where x = Oand 0 S B S z. By a well known 
property of uniform magnetic shells which remain valid in any number of dimensions, 
the stresses r,, and 7», remain continuous across the disc. However, they become infinite 
in the neighborhood of its rim like (y — b)[a” + (6 — y)*]’ and a[a° + (6 — y) 
respectively. The connection of the formulae (6.1) and (6.2) with elliptic integrals 
analogous to those given in [4] and [5] will be given elsewhere. Further applications of 
‘““seneralized electrostatics” are given in [10]. 

The author wishes to thank Dr. H. F. Weinberger for valuable discussion on the 


material of this paper 
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INDUCED MASS WITH FREE BOUNDARIES* 
sy GARRETT BIRKHOFF (Harvard University) 


The concept of induced mass, and some of its properties, are extended to the case 
of an incompressible liquid having a free surface. The usual’ treatment of the case of a 
non-viscous fluid extending to infinity must be considerably changed to do this. 

1. Minimum principle. Let F be a region filled with an incompressible liquid, bounded 
in part by a wetted wall W, and in part by a free surface S at constant pressure p, . The 
region R moves with the liquid, and may extend to infinity in some directions. 

Suppose the fluid accelerated from rest, by an acceleration of W whose inward normal 


component toward the liquid is an arbitrary function f(x) of position. Letting u = 
u(x; ¢) denote liquid velocity and a = du/dt denote acceleration, clearly 


Div a = 2d(du,/dl)/dx, = A(Tdu,/dx,)/dt = 0, (1) 


by incompressibility. Similarly, the Navier-Stokes equations with gravity neglected are’ 
: ] : i | { y g 


Du;/Dt = —dp/pdx; + vT>u; « (2) 
Since u = 0 initially, the znztial acceleration therefore satisfies a; = Du,/Dt = —dp/pdz; , 
or, setting A = (po — p)/p, 
a= VA, where Pp = Po — pA, initially. (3) 
Combining (3) with (1), we get 
V°A =0 in R. (4) 
The free surface condition is simply 
A=0 on S. (5) 
Finally, by continuity, we have 
Gnormat = OA/On = f(x) on WW. (6) 


*Received July 16, 1951. 
1Given in [1], Chap. V; in [3]; and in Chap. VI of Lamb’s Hydrodynamics. 
2Here and below D/Dt is the substantial derivative 0/dt + ux 0/dz; , while 0/8 n denotes the inward 


normal derivative on the surface of the liquid. 
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By potential theory,’ conditions (4) to (6) uniquely determine A; we shall call A the 
acceleration potential associated with the given acceleration of W; by (3), A determines p. 


THEOREM 1. The acceleration kinetic energy 


T => [/| VA:VAdR (7) 
J 
is minimized by the free surface condition (5), relative to all other volume-conserving 
flows in R satisfying (6). 
Remark. The acceleration kinetic energy is half the second time derivative of the 
ordinary kinetic energy. 
Proof. Let VA + b be any other volume-conserving flow satisfying (6). Then 


Div b = Divb + V’°A = Div (VA + b) = 0 


Consider now the expanded acceleration kinetic energy 


ny 


T=s5/ (VA+b)(VA +b) aR 


oo “R 


iar lg l j 
=T,+p/ (VA-b)aR+i0/| (b-b) aR. 
YR =“ %/R 
Since the last term is positive unless b vanishes identically, the theorem will be proved 
if we can show that the middle integral is zero. But since Div b = 0, clearly 


Div (A-b) = A Div b + (VA)-b = VA-D. 


Hence, by the Divergence Theorem, letting b,, denote the outward normal component of b, 


. sy 


| (VA-b)dR = —| Ab,dS— | Ab, dS. (8) 
J Jn Js 
Since VA + D satisfies (6) on W, b, = O on W and the first surface integral in (8) 
vanishes. By (5), the second surface integral vanishes; hence the proof is complete. 

When S is void, our result reduces to the classical case ({2], p. 84). There is however 
no relation between the kinetic energy for sfeady flow with a free boundary, which is 
always infinite in an infinite non-viscous stream, and virtual mass. In this respect, the 
case of flows with free boundaries is unlike the classical case. 

2. General interpretation. The acceleration potential A described in Sec. 1 has a 
very general interpretation. Let a fluid motion, defined by an “undisturbed” velocity 
field U = U(z; #), be altered at time ¢ = 0 by an additional instantaneous normal ac- 
celeration f(x) of W. Because of incompressibility, we have for the perturbed velocity 


field U = u, 
Div u = 0, whence Div a = 0, (9) 


where a = 0u/di denotes the additional acceleration. 

3For existence and uniqueness theorems, see O. D. Kellogg, Potential theory, pp. 218, 315. Although 
such theorems have been proved only for a restricted class of infinite regions, there is no reason to doubt 
their general validity. 
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Subtraction of the Navier-Stokes equations for U from those for U + u gives 


Ou; ou, , OU; Ou; 1 Op 
= —U, _ _ em Sm aan 


= = _ ree 
at Ox, * Ox, ' Ox). p Ox; 


+ vV*u; . (10) 


Here p denotes the increment in the pressure field due to the additional acceleration 
of W; moreover, this is true even if gravity is considered. Since the perturbation begins 
at ¢ = 0, clearly u(x; 0) = 0, whence (10) reduces to 


a; = —dp/pdz,; , at t= 0. (11) 


Since the rest of the reasoning leading to (4) to (7) applies, we get the following result. 
THEOREM 2. The instantaneous pressure distribution required to accelerate a moving 
incompressible fluid is the same as if the fluid were at rest. 
The preceding result can be extended to the case of an “impulsive” velocity change, 
by integrating (10) over a short interval of time, which is then allowed to tend to zero. 
It is readily seen that if u is uniformly bounded and ff V~u dR dt = o(1) under these 


circumstances, and if the normal impulse per unit area is defined as p* = Lim,, fo‘ pdt, 
we get the limiting analog of equation (3), 
u; = —dp*/pdx, on W. (3*) 

From (9), the fact that p* = 0 on the free surface S, and continuity we get similarly 
equations analogous to (4) to (6). From these we deduce the following result. 

fHeoreM 3. The impulse p*(x) per unit area required to produce an additional normal 
“impulsive velocity” u, = f(x) of W, is —pA, where A is the acceleration potential 
of Sec. | 

In the case of the irrotational motion of a non-viscous fluid, since V*u = 0, the 
assumption about Y~u (which is hard to prove rigorously) is superfluous. 


3. Induced mass tensor. To define an induced mass tensor, lei W = W’ + W” con- 
sist in part of the wetted area W’ of an accelerated missile, and otherwise of the fixed 
walls W” of a container. The case of a ball, floating on the surface of the water in a pail 


(Fig. 1) is typical. 














We let A’, A’, A® denote the acceleration potentials in R for unit translations of W’ 
parallel to the axes; let A*, A°, A° denote those for unit rigid rotations of W’ about the 
coordinate axes; in all cases we assume 0A/dn = 0 on W”, and A = 0 on S. We then 
define the symmetric 6 X 6 induced mass tensor (matrix) || 7’, ||, as in the ordinary 


case ({1], p. 154), by 


Tu = p Il VA'-VA' dR = Ty (12) 








84 NOTES [Vol. X, No. 1 


Our aim is to show that the 7',, have most of their familiar properties. Although the 
preceding definition of the 7, was given in principle by L. I. Sedov [4], he did not 
deduce the properties proved below. 

First, note that the diagonal components 7',, of induced mass, satisfy the conditions 
of Theorem 1. From this fact the following result follows immediately. 

Coro.uary 1. The diagonal components 7',, of induced mass are increased if either 
(i) the region R occupied by liquid is increased (i.e., the free boundaries are pushed out), 
or (ii) free surface area is replaced by container walls or wetted missile area. 

Corouuary 2. Let a volume AV of liquid be replaced by missile. The new translation 
induced mass 7%, satisfies the inequality. 

Tin = Ta — pAV (13) 

Proof. Let R’ = R AV denote the reduced volume occupied by liquid. Consider 
the acceleration field of R, under which AV is accelerated as a rigid body under h- 
translation, while R’ is given the acceleration corresponding to Tj, . By Theorem 1, 
Ti, + pAV will exceed 7, , proving formula (13). This result can be extended to rota- 
tional components, if pAV is replaced by the appropriate rigid moment of inertia. 


Again, ve have by Green’s second identity 
vas f 0A” , F , OA A i OA’ : 
—7 =e ft 2 > dS + p | A — dS + p | A’ — das 
j on Jn on Jw on 
By (5), A* = 0 on S, and so the first summand is zero; by (6), dA"/dn = 0 on W”", 
hence the last summand is zero. 
Using 3) and reversing signs, we therefore have 
i OA ; 
1 = | Dp _ Po) . ds, (14) 
™ on . 


“W 


where p, is the free surface pressure, and p, is the pressure on W’ under unit h-accelera- 


tion from rest. Also on IW’, dA‘/dn dx*/dn [k = 1, 2, 3] and 8A*/dn = 2,02;/dx, — 
2,02,/0x, ete., by definition. Hence, cancelling out the effect of the constant pressure 
Py over the missile, we have 


TureoreM 4, The tensor component 7’, represents the total h-component of initial 
pressure fo "CE required to produce a unit k-acceleration of the missile. 

In view of the linearity of (3) to (6) in A and p, we deduce immediately the following 
corollary 

Coro.uary 1. If, at any instant, the configuration and flow field are given, the in- 
stantaneous effect of an additional acceleration of the missile, with components a, , --- 6, 
is to produce components } , T',,a, of reaction in the liquid. 

This result will be compared with experimental data in See. 5. 

4. Momentum interpretation. If there is no container, so that W = W’ and W”’ = 0, 
momentum interpretations of the 7’,, are also possible. I shall show this directly from 
momentum considerations. (To make this rigorous, one must consider, as in [1] or [3], 
the convergence of the momentum integrals involved.) 

Let the missile, supposed of finite diameter, be given initial acceleration from rest, 
and let C by any cylinder parallel to the 2z,-axis, whose finite cross-section contains the 
wetted area W = W’ of the missile. The total h-momentum* in C is finite; moreover, 





4By h-momentum, we mean linear momentum parallel to the x-axis. 
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since the liquid is at rest initially, the rate of convection of fluid h-momentum out of 
C is zero. Since the fluid is non-viscous, the shear stresses across this boundary are also 
zero; hence the rate of total transfer of h-momentum across the boundary of C is zero. 
We conclude that the time rate of increase of h-momentum of the liquid in C, is equal 
to the h-component of pressure force. 

A corresponding result holds for moments of momentum; in this case C must be a 
solid of revolution containing W, with the axis about which the moment is taken, for 
axis of symmetry. Combining with Theorem 2, we get the following result. (We define 
an “h-curve” to be a straight line parallel to the x,-axis, if h = 1, 2, 3; and to be a circle 
perpendicular to the x,_3-axis, with center on the x,_;-axis, if h = 4, 5, 6.) 

TuHeoreM 5. Let a missile be given a unit k-acceleration from rest, in a liquid bounded 
by the missile wetted area and a free surface. Then the rate of increase of the h-com- 
ponent of liquid momentum in any region C bounded by h-curves, which contains WV, 
is exactly Th initially. 

5. Application. we can apply the corollary of Theorem 4 to the case of a missile B 
travelling vertically in a cavity, as in Fig. 2. The reaction to acceleration of B will be 
at least as great as if the liquid were confined to the underside of a plane bounding the 


wetted area of B (indicated as a dashed line in Fig. 2). 





In this case, as first noted by von Karman [6], the acceleration potential can be 
obtained by symmetry, using the reflection principle. It follows that the instantaneous 
inertia opposed by the liquid to the acceleration of B, should be between 50 per cent 
and 100 per cent of that offe red if there were no cavity, and the same area were wetted. 

This is not necessarily the same as the change in the cavity drag coefficient due to 
steady acceleration or deceleration (ef. See. 1, end). In fact, data obtained at the Naval 
Ordnance Laboratory’, indicate that this change probably corresponds to less than 25 
per cent of that occurring if there were no cavity. 

The preceding results also apply to the case of impact on water, for which the mo- 
mentum interpretation of §4 may be of interest. In this case, a rigorous lower bound to 
the loss of energy at impact is obtained by Theorem 1. Unfortunately, about 50 per 
cent of the energy lost at impact is presumably absorbed by the energy of compression 
(and radiated as acoustic energy); hence this bound is excessively low. 

However, as the impact phase has been extensively discussed by other authors 
({4], [5], [6]), especially in connection with the landing of seaplane floats, we shall not 
discuss it further here. 

5A. May and J. Woodhull, ‘The virtual mass of a sphere entering water vertically”, Journal of 
Applied Physics 21 1285-9 (1950). 
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6. Correction for hydrostatic force. Since gravity has been neglected above, it is 
interesting to have a rough estimate of the effect of gravity on the pressure exerted by 
a liquid on a missile moving through it with wetted area W’ trailed by a cavity (Fig. 2). 
We suppose the liquid incompressible, and bounded by W’, container walls W’’, and a 
free surface S. The additional instantaneous acceleration b due to a vertical gravity 
field with intensity g satisfies b = gVB, where V*B = 0, B = y (depth coordinate) 
on S, and dB/dn = 0on W’ + W” = W;; the associated hydrostatic pressure is pg(y — B). 

For given boundary configurations S and JW, the resulting “hydrostatic acceleration 
potential’’ gB can be most easily estimated using an electrolytic tank, and the results 
interpreted in terms of the dimensionless cavity buoyancy coefficient 


hydrostatic bouyancy force 





Cg @ eee : ie Rae Pm 
pg X mean depth X horizontal projection of W 


In this way, Cy was estimated’ for three two-dimensional cavity flows, having profiles 


similar to that of Fig. 2. The cavity buoyancy coefficients averaged about 25 per cent. 
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ON THE NON-UNIQUENESS OF PERIODIC SOLUTIONS FOR AN 
ASYMMETRIC LIENARD EQUATION* 


By G. F. D. DUFF anp N. LEVINSON (Massachusetts Institute of Technology) 
The following result has been stated by H. Serbin [5, Theorem II]. Let f(x), g(x) 


be continuous for —7 < x <o, and let 


fiz) <0, -xi<r<n, (1.0) 


f(x) > 0, <2, ae (1.1) 


, 


where x; > O and z, > O. Let 


1 


| f(x) dx > 0 (1.2) 


“0 


xg(x) > 0, r#0 (1.3) 


*Received Aug. 20, 1951. This paper was written in the course of research sponsored in part by 


the Office of Naval Research. 
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a) 
0 


f(x) dx or f¢ g(x) dx diverge. Then the equation 


x+ f(x)x+ g(t) =0, x = dx/dt, (1.4) 


and let either | 


. . . . ° 
has a unique periodic solution. 


In case g(x) > and f(x) is even (1.4) has been treated by Liénard [3]. In case 
f(x) is even and g(x) is an odd function (1.4) has been treated by Levinson and Smith 
[2]. Here an example will be given with g(z) = x and with f(x) meeting the requirements 
(1.0), (1.1) and (1.2) where at least three periodic solutions for (1.4) exist thus showing 
that the result of Serbin is incorrect as stated. 

We note that if f(x) satisfies the required conditions so does ef(x) for any « > 0. 
The equation to be considered here is 





rt+ef(sxr+x2=0 (1.5) 
where f(z) is a polynomial. If 2 = r cos @and x = —r sin 6 then (1.5) yields r = er 
f(r cos @) sin® 6, @ = 1+ ef(r cos @) sin 6 cos @ from which follows 
dr rf(r cos @) sin’ 6 
© i oe. (1.6) 
dé 1 + ef(r cos @) sin @ cos 6 

For each solution of (1.6) r = r(@) periodic of period 27 there corresponds a periodic 


solution of (1.5). 

Since with « 0 the solutions of (1.6) r = const., are bounded and since the right 
side of (1.6) is analytic in r, 6 and e¢ it follows from the successive approximations existence 
theorem, as indicated for example in Lefschetz [1, Chap. II] that the solutions of (1.6) 
with r = p for 6 = O are given by r = H(8@, p, «) where H is analytic in (6, p, €) for any 
finite range such as 0 < @ < 2x, 4 < p < 5, and | €| sufficiently small. Thus H can be 


expanded in powers of ¢ yielding 


r = H,(0, p) + eH,(0, p) + ey (8, p, €) 


where | y| is uniformly bounded for the above range of 6 and p for | € | small. Using 
r = H in (1.6) it follows, equating powers of ¢ on each side of the equation, that 1,(@, p) = 
p and 
H6, p) = | pf(p cos a) sin’ o de. 
It 
F(p) = | pf(p cos a) sin’ o do (1.7) 
then 
H(2r, p, €&) — p = €F(p) + &y(2r, p, ©). (1.8) 
If H(2z, p, €) — p vanishes for some p then as already stated (1.5) has a periodic solution 
passing through the point x = p, x = 0 in the phase plane of the variables (2, z ). 


Let po > O and let F(p,) = O and let F(p) change sign at p = p, . Then for a given 
small 6 > 0, it follows from (1.8) that H(2z, p, «) — p has opposite signs for p = po — 6 
and p = p, + 64 providing « is small enough. Since H(2z, p, «) — p is continuous in p 
this implies that 7(2z, p, «) — p vanishes for some p between py — 6 and po + 6. 

It follows therefore that if it is shown that F(p) in (1.7) has 3 positive zeros and 
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i? 2) 
wa 


changes sign at each zero then for small ¢ (1.5) has at least three periodic solutions as 
stated. [4] 
Let 
f(a) = Agu® — Ana* + Aix? — Ay — Cx (1.9) 


where A; and C are constants. Let 


I, = | sin’s cos o do, k= 0; 2.2: S. 


“0 


Then J, > 0 and from (1.7) it follows that 
F(p) = AsIz3p’ — Aslop’ + Aylip* — Aolop. 
Let Al, = 1, Aal. 14, A,J, = 49 and A,J, = 36. Then all the A; are positive and 


F(p) = p(p’ — 1)(p — 4)(e° — 9). 
Thus F(p) vanishes at p = 1, 2, and 3 and changes sign at each of these zeros which 
proves that (1.5) has at least three periodic solutions for small e. 

Since (1.2) is assured by A; > 0 in (1.9) it remains only to show that f(x) satisfies - 
(1.0) and (1.1). Since the A; are already fixed and since A; > 0 and A, > 0 it follows 
easily that if C is large enough then Cr = A,x° — A.x* + Ax” — A, has only one 
negative root, which is near z = 0 and will be designated by —2{ , and only one positive 
root, which is large and will be designated by x, . Clearly f(x) now satisfies (1.0) and 
(1.1) and the example is complete. 

With polynomials of higher degree the existence of more periodic solutions can be 
established. (The argument of Serbin fails in that he does not recognize the restricted 
range of the monotonicity of ¢, in Lefschetz’s notation, as stated by Lefschetz [1, p. 193] 
and to which Serbin refers.) 
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THE COMPRESSIBLE FLOW CORRESPONDING TO A LINE DOUBLET* 
By J. W. CRAGGS (St. John’s College, Cambridge) 

Introduction. The use of the hodograph equations for irrotational compressible flow 
in two dimensions leads, as noticed by Chaplygin’ and others, to a method of con- 
structing compressible flow patterns by reference to similar problems in incompressible 
flow. For certain elementary flows, as for example the flow due to a line source, a line 


*Received December 11, 1950. 
1Chaplygin, On gaseous jets, Ann. Sci. Moscow Imperial University, Math.-Phys. section 21, (1904). 
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vortex or a line doublet, the method is particularly simple. Ringleb* has investigated 
these cases, but his diagrams for the doublet do not illustrate what is perhaps one of 
the most interesting things about the corresponding compressible flow, the type of limit 
lines and singularities which occur. In the present paper special attention is paid to 
these points. 
The solution in the hodograph plane. For incompressible flow, the complex potential 
for a line doublet is given by 
w= —1/z, ge *’ = dw/dz = 1/2’, 
where q is the magnitude and @ the inclination to the real axis of the velocity. 


Eliminating z we obtain 


whence y = q'” sin 0/2, (1) 


where y is the stream function. The two branches of (1) correspond to the two half- 

planes, on each side of the doublet axis. If only the upper half-plane is considered the 

hodograph plane is once covered, but is cut from the origin along 6 = 0 to q = 1. 
That solution of the hodograph equations of compressible flow which corresponds 


to (1) is 

v= q' f(r) sin 6/2, (2) 
where T= QQ 
Q being the maximum velocity, corresponding to zero density, and 

f(r) = F(a, b; 3/2; 7) (3) 
where F is a hypergeometric function with 

ab = —3 8, a + b= (1 = 28) 2, 
being the reciprocal of y — 1, where y is the adiabatic index. 
The hypergeometric series (3) converges rapidly and is easily calculated for values 

of 7 < 1. Nearr = 1 the expansion in terms of 1 — 7 may be used.* Thus the flow pattern 


in the hodograph plane is found. 
Reversion to the physical plane. The coordinates in the physical plane, x, y, are 
given by 
dx = (d¢/q) cos 8 — (po/p)(dy/q) sin 
(4) 


and 


(dp/q) sin 8 + (p0/p)(dy/q) cos 8, 


dy 
where 
q(d¢ /Aq) - —(po/p)(1 _ M’)(dy/d8), 
(5) 


(06/00) = (po/p)q(d¥/0q), 
*Ringleb, Exact solutions of the differential equations of an adiabatic gas flow. Z.A.M.M. (1940). 
*The form of the alternative expansion may be found in any standard text on hypergeometric 


functions. 
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M being the local Mach number, given by 
M? = 26r/(1 — 7) 


and (p/p) the density ratio, given by 
2 —28 
(po/p) = (l— 1) ™. 
Equations (2), (4) and (5) lead to expressions for the derivatives of x and y with 
respect to q and @. Integration of those expressions then gives 


a(pq'’’/po) = 27(df/dr)(cos 36 — 4 cos 30) + f cos 40 
and (6) 
y(pq'’*/po) = 2r(df/dr)(sin 46 — 4sin 36) + fsin 30. 


The flow in the hodograph plane may be found by plotting streamlines y = const. 
The streamlines in the physical plane are then obtained by substituting corresponding 
values of g and @ in (6). 

Singularities of the transformation. Now we investigate the singularities of the trans- 
formation. The transformation from the hodograph to the physical plane ceases to be 
regular when 

J = d(q, 6)/d(a, y) 
vanishes or becomes infinite.* The Jacobian reduces to 


J = —(p/po)'9/{a ve + (1 — M)* yi}. 


Zeros of J can occur only when g = 0 or y, or ¥y ~~. Now f and df/dr are always 
finite so the last two possibilities are eliminated by (2). As q — 0, f —~ 1 and J — 0; 
this corresponds to a singular point (the point at infinity) in the flow plane, and a branch 
point in the hodograph plane. The singularity is due to the symmetry of the flow about 
the x and y axes, which implies the same velocity vector at (x, y) and (—xz, —y). This 
has already been dealt with by the cut in the hodograph plane. 

The singularity corresponding to J —© is of more interest. J © implies gy, + 
(1 — M’)y¥; = 0. This can only be satisfied for supersonic velocity, M > 1. Put a’ 


M*~ = 1, then the condition is 
QW, + av, = 0, (7 
that is 
(2rf’ + 4f)q'”? sin 46 + 3aq'f cos 46 = 0, 
whence 
a cot 46 (4rf’/f) + 1. (8) 


The singular line in the hodograph plane, given by (7), is a double loop curve touching 
the sonic circle, M = 1, at 6 = 0, crossing itself at 6 = 180°, 4r(f’/f) + 1 = 0, (r = 
.451, approximately) and touching the limit circle, g = Q at 6 = 180°. It is shown in 
Fig. 1 together with plots of typical streamlines in the hodograph plane. 


7c) s, The breakdown of the hodograph transformation for compressible flow in two dimensions 


Liveypy 


Camb. Phil. Trans. 44, 1948. 
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The limit line in the flow plane may now be found by substitution in (6) of (8). 
Corresponding to the double point there are zeros of dx/dr and dy/dr taken along the 








q=Q 


Fia. 1. Hodograph plane. 


curve, and each limit line is cusped. It is easily deduced from (2) and (5) that at this 


point 

Vo = Vo = b, = bo = Woo = bee = Goo = O 
whilst 96.0 = (po/plaWoo 
and TVoa = a Woo , 


thus the ratio d¢/dy corresponding to given dq/d@ at this point is 


dp _ 2a°(p./p)q(dq/dé) 
dy a (dq/d@)’ + ¢ ° 


But on a limit line gd@/dqg = +a, whence d¢/dy = +a(po/p). It follows from (4) that 


] os 6 = sin 8 
2 a = For ft (Oe 88. 
dx asin + cos @ 


Next, since yg, = 0 and y,, , Wee are of the same sign, we see that the point is an ex- 
tremum of ¥ and hence that lines Y = const. encircle the point and no streamline passes 
through it. The singularity is of order two (see ref. 3). 

The flow in the physical plane. We may now see the main characteristics of the flow in 
the physical plane. First we draw in the limit lines, shown chain dotted in Fig. 2. Each 
of the four branches leaving the point B is the join of two sheets of a Riemann surface, 
the four sheets of which correspond to the regions a, 8, y, 6 of Fig. 1. The regions 8 
and 6 are separated by the cut @ = 0. 

The region a corresponds to a sheet covering the whole of the diagram above the 
lower pair of limit lines, BP, BQ. Along BQ it joins a sheet (8) bounded by BQ and BR. 
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The sheet y extends from BR to BS and the fourth, (6) bounded by BS and BP, is 
symmetrical to the second, and joins the first (2) along BP. A streamline, like that for 
y = 0.92, which is closed in the hodograph plane, passes over the four sheets in order. 





—-—-—LIMIT LINE 


‘ ————STREAMLINE SHEET OL \ 
----- OTHER SHEETS \ 
j \ 
| j 
Fic. 2. Physical plane. 
The stream line for Y = 0.5, which is also shown, never reaches the sheet y but goes 


off, after reflexion at each point where it meets the limit lines, to infinity with speed at 
infinity of Q. 


TaBLeE I. The Limit Line. 


T 6 zx y 
.16 0.00 1.278 .000 
18 58.36 1.032 .684 
.20 87.48 .715 .906 
22 107 .00 .555 . 986 
.24 121.00 .422 1.015 
26 132.12 .312 1.020 
.28 141.24 . 230 1.014 
30 148.48 . 167 1.007 
.35 162.42 .062 .997 

40 72.36 .014 .959 
.45 179.20 .001 .948 
.50 184.14 .010 .955 
.55 187.42 .041 .975 
.60 190.02 .095 1.003 
.65 190.42 . 184 1.063 
.70 191.10 .320 1.170 
.75 192.00 589 1.314 
.80 190.44 .965 1.600 
.85 189.56 1.859 2.060 
.90 188.02 4.280 3.760 
.95 183.44 18.3 15.4 

1.00 180.00 +o +0 
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TaBLE II. Streamlines 


y¥ = 0.5 yy = 0.9 ¥ = 0.92 

T x y zx y zx y 
.00 

.02 . 769 1.618 

.03 .92 1.347 

.04 .962 1.190 

.05 .975 1.087 

.10 .961 .860 

12 .952 .824 

.14 .945 .801 

.16 .942 . 788 

.18 .940 .78 

.20 .940 . 784 

22 .943 .790 

24 .947 .801 .065 1.017 

.26 .954 .816 .105 1.007 

.28 .962 . 836 .123 1.001 

.30 .973 .859 .131 .997 .042 . 983 
.35 1.009 .938 .124 1.000 .068 .976 
.40 1.063 1.051 .096 1.019 .054 . 983 
.45 1.138 1.207 .049 1.049 .019 .998 
50 1.241 1.424 — .018 1.091 — .031 1.018 
55 1.387 1.731 — .108 1.148 — .904 1.042 
60 1.599 2.181 — .230 1.219 — .163 1.067 
65 1.922 2.822 — .397 1.318 — .234 1.092 
70 2.439 3.956 — .627 1.441 — .126 1.091 
75 3.333 5.986 — .959 1.595 

.80 5.09 9.97 — 1.426 1.798 

.85 9.25 19.59 —1.766 2.022 

.90 22.8 51.85 

.95 117.4 282.6 

1.00 


For clarity only y = 0.5 and y = 0.92 are shown in the diagrams, 


BOOK REVIEWS 


Mécanique des milieux continus et déformables. By Maurice Roy. Gauthier-Villars, Paris, 
1950. Volume I, xxii + 363 pp., 2800 fr. Volume II, xii + 350 pp., 2300 fr. 


M. Roy presents an account of the mechanics of fluids and elastic solids that is based upon the 
application of thermodynamic principles to continua. The work covers substantially the content of a 
second year course of lectures at I’Ecole Polytechnique, and the object is to give a synthetic account 
of the fundamental principles and methods of this mechanics bearing in mind applications to engineering. 

Each volume comprises two parts together with extensive appendices; the four parts are entitled 
“thermodynamics and mechanics of continuous and deformable media”, ‘theory of elasticity’’, “equilib- 
rium and flow of fluids’, and “theory of machines’. The principal contents of these parts are described 
briefly as follows. Part I. “Definitions and fundamentals of deformations”; finite and infinitesimal 
strains, Eulerian and Lagrangian co-ordinates, acceleration, equation of continuity. “Fundamental 
thermodynamic principles and equations’’; thermodynamic variables, first and second laws of thermo- 
dynamics, interna] potential, general thermodynamic equation. “General theorems of mechanics’’; 
stability of equilibrium, reciprocal theorems. “Applications to deformable media and elastic solids’’; 
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distinction between elastic solids and fluids, stress, infinitesimal deformation of elastic solids, Hooke’s 
law, coefficients of elasticity, equations of equilibrium, isotropic and homogeneous solids, elastic con- 
stants. ‘‘Application to viscous and non-viscous fluids”; thermodynamic definition of a fluid, viscosity, 
Navier’s hypothesis, dissipation function, equations of continuity and equilibrium and state, Reynolds 
and Mach numbers. Part II. ‘‘Generalities, problems and equations and compatibility”; elastic equilib- 
rium, Clapeyron’s formula, Castigliano’s theorem, elastic limit. “Plane elasticity”; plane stress and 
strain, Airy’s stress function, use of complex variables, examples. ‘“Three-dimensional elasticity”’ ; simple 
tension problems, simple flexure of beams, pure torsion, Saint-Venant’s problem, beam theory, theorem 
of three moments. ‘Thin bodies and elastic bars”; elastic properties of bars, stability. ‘Small motions 
of three-dimensional elastic bodies”; dilatational and rotational waves. Vibration of thin bars is con- 
sidered in the appendix. Part III. ‘‘Statics of fluids’; equilibrium conditions, stability of floating bodies. 
“Dynamics of non-viscous fluids’; Eulerian equations of motion, Helmholtz’ circulation theorem, 
Torricelli’s and Saint-Venant’s formulae, d’Alembert’s paradox. ‘Plane irrotational flow of incom- 
pressible fluids’; velocity potential, stream function, complex potential, singularities in flow field, 
Kutta-Joukowsky hypothesis, Blasius’ formulae, Kutta-Joukowsky theorem, aerofoil theory, images, 
dead-water regions and cavitation, jets, vortexes, Karman vortex streets. ‘“Three-dimensional irro- 
tational flow of incompressible fluids”; sources and sinks, images, Rankine solids. ‘““Three-dimensional 
rotational flow of incompressible fluids’; free and bound vortex sheets, velocity induced by vortex ring, 
Prandtl’s lifting-line theory, air-screw theory. ‘‘Wave propagation’”’; propagation of small disturbances 
in compressible fluid, shock waves, Hugoniot relations, stationary shock, waves on a canal. ‘“‘Com- 
pressible fluid flow’’; potential flow, one-dimensional motion and blast waves, rectilinear singularities 
in flow field, plane subsonic potential flow with the use of hodograph plane, Prandtl-Meyer expansions, 
characteristics, simple waves. ‘‘Dynamics of viscous fluids’; slow motions, laminar motions, turbulent 
motions, hydraulics, stability of laminar flow, Prandtl’s boundary layer theory, transition, separation. 
Part IV. ‘‘Case of a periodic regime’’; fundamental theory, application to turbo-machinery. ‘‘Case of a 


varying regime”’; application to hydraulic turbines, stability of regime. 

The reviewer commends the author for the clear and stimulating treatment of his subject. The 
work is a useful contribution to the literature on the mechanics of continua. The reviewer regrets that 
the tensor notation, with the summation convention, is not used; as Harold Jeffreys has remarked, this 
notation, although not a necessity, is a great convenience and, when cartesian axes are employed, is not 
complicated. If discussion were given of this notation, then, in the reviewer’s opinion, the presentation 
of the fundamental analysis would be considerably improved. The list of contents, although very de- 
tailed, does not avoid the need for an index. A rather surprising feature of the work is that comparatively 
few references are given. 


H. G. Hopkins 


Fourier transforms. By Ian N. Sneddon. First Edition. McGraw-Hill Book Company, 
Inc., New York, Toronto, London, 1951. xii + 542 pp. $10.00. 


The aim (stated in the preface) of this book is to present the theory of Fourier transforms and 
related topics in a form suitable for the use of students and research workers interested in the boundary 
value problems of physics and engineering. The book is written for those interested in the application 
of the theory rather than the theory itself. 

The theory of Fourier transforms, of Hankel transforms and of finite Fourier and finite Hankel 


transforms is covered in the first three chapters (90 pages) of the book. In this section of the book there 
appears, among other things, discussion of integral transforms in general, Fourier integra] theorem, 
LaPlace transform, Mellin transform and multiple Fourier transforms, Hankel inversion theorem, 
Hankel transforms of the derivatives of a function, relation between Hankel and Fourier transforms, 
dual integral equations and the Mellin transform. The Wiener Hopf technique and the delta function 
are discussed and used in later chapters. The transform theory discussed in the first fifth of the book is 
applied to physical problems of sufficient difficulty to be of value to the research worker. 
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The chapter headings in that part of the book devoted to applications are as follows: 4: Theory of 
Vibrations, 5: The Conduction of Heat in Solids, 6: The Slowing Down of Neutrons in Matter, 7: 
Hydrodynamic Problems, 8: Applications to Atomic and Nuclear Physics, 9: Two Dimensional Stress 
Systems, 10: Axially Symmetric Stress Systems, Appendix A; Some Properties of Bessel Functions, 
Appendix B: Approximate Methods of Calculating Integral Transforms, Appendix C: Tables of Integral 
Transforms (Fourier, LaPlace, Mellin, Hankel, Finite Fourier, Finite Hankel). 

The text has many virtues, one of which is the development of the discussion in each chapter from 
fundamentals; something that students will certainly appreciate. The book seems to be a very valuable 
contribution to the literature of physics and applied mathematics. 

Roun TRUELL 


Introduction to the study of aircraft vibration and flutter. By R. H. Scanlan and R. 
tosenbaum. The Macmillan Company, New York, 1951. x + 428 pp. $7.50. 


This text is an introductory treatment of questions of mathematics, structural theory and aero- 
dynamics which are of importance for the analysis of aircraft vibrations, forced or self-induced. It is 
written for senior or graduate engineering students but should also be useful to others wishing to learn 


the rudiments of the problem of “flutter’’. 
E. REISSNER 


Quantum mechanics of particles and wave fields. By Arthur March. John Wiley & Sons, 
Inc., New York, and Chapman & Hall, Limited, London, 1951. x + 292 pp. $5.50. 


It might seem difficult to justify another book on quantum mechanics, yet this book by Arthur 
March seems to cover concisely, and well, quite a few topics normally omitted in the elementary texts 
on quantum mechanics. The chief virtue of the book seems to be that it helps to bridge the gap between 
elementary texts and the current literature. It will be found that about half of the book deals with 
relativistic quantum mechanics. The book is concerned primarily with the exposition of general principles 
and it omits methods of calculation and many applications. The result seems to the reviewer to justify 
the omissions. 

The last chapter deals with the idea of universal length, an idea apparently first set forth by the 
author. While the ideas of fundamental length do not seem to be a part of current theory, calculations 
based on these ideas do, in some cases at least, give agreement with experiment. 

Roun TRUELL 


Dimensional analysis and theory of models. By Henry L. Langhaar. John Wiley and Sons, 
Inc., New York, and Chapman & Hall, Limited, London, 1951. xi + 166 pp. $4.00. 


Dimensional Analysis and the Theory of Models are subjects which have been in common engi- 
neering use for many years. The basic ideas are so elementary that adequate practical results can 
frequently be obtained without a very rigorous understanding of the basic ideas on which the process 
rests. In this book, the author attempts to bring together a complete picture of these subjects. The 
first two chapters in fact, are devoted to a discussion of dimensional units, the elementary development 
of Buckingham’s Theorem and a series of applications to elementary problems. The third chapter is 
devoted to the systematic procedures by which dimensionless products may be deduced. Thus, the 
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first three chapters constitute what is sometimes regarded as the whole of dimensional analysis based 
upon intuitive arguments. 

In chapter 4, Algebraic Theory of Dimensiona] Analysis, the author has brought together a simple 
but concise treatment of what is now known of the mathematical foundations of dimensional analysis. 
In this chapter, Buckingham’s Theorem is proved without appeal to the possible expansibility of various 
functions into power series. After thus providing the sound foundation for dimensional analysis, the 
book takes up in succession a series of specific examples. The items considered include Ship Propeller, 
Underwater Explosions, Gas Flow, Archery Bow, Bending of Materials, Small Deflection of Elastic 
Structures, Free Vibrations, Turbulence, Boundary Layers, Run off From a Water Shed, Centrifuga- 
Compressors, Heat Transfer in a Pipe, Condensation, Transient Heat Transfer, Thermistors, Piezo- 
electric Receivers, etc. Before taking up this last problem, a careful treatment is given to electrical and 
magnetic units from the point of view of dimensional analysis. 

The treatment of model theory from the point of view of model—prototype similarity starting 
with the differential equations is considered in chapter 10. 

The book makes no pretense to contribute to knowledge on the way which a person must proceed 
to select the appropriate variables and appropriate basic dimensions for a given problem. To be sure, 
the numerous examples would assist the student in attaining an understanding of how to go about it, 
but the basic problem still remains. Thus, no clear understanding is yet available to guide the un- 
initiated with this task. The reviewer does not consider the problem satisfactorily disposed of by the 
statement that experience is the only guide. This is, of course, true at the present time but is to be 
regarded as an unsatisfactory temporary situation. 

For anyone uninitiated in the ways of dimensional analysis and the use of models or anyone wishing 


an excellent systematic and thorough treatment of the subject, this book is excellent. 


Howarp W. Emmons 


Tensor Calculus. By J. L. Synge and A. Schild. University of Toronto Press, Toronto, 
1949. xi + 324 pp. $6.00. 


This book presents tensor calculus in a manner well suited both to those interested in an abstract 
study of geometry of spaces with general dimensionality, and to the increasing number of those in- 
interested in the applications of tensor calculus to mathematical physics and engineering. This book is 
an excellent classroom text, since it is clearly written, contains numerous problems and exercises, and 
at the end of each chapter has a summary of the significant results of the chapter. 

Chapter 1 deals with the definition of a space, and the definition and algebra of tensors. Chapters 
2, 3 and 4 deal with Riemannian spaces of general dimensionality whose fundamental quadratic forms 
are not necessarily positive definite, and include the concepts of distance, angle, geodesics, tensor 
derivatives and the Frenet formulas. There are also extensive treatments of curvature of space and 
cartesian tensors. 

Chapter 5 contains a study in terms of tensors of the dynamics of a single particle and of a rigid 
body. There is also here a discussion of the geometry and topology of various configuration spaces of 
dynamical systems. In Chapter 6 there appear concise developments of the fundamental equations of 
the mathematical theories of electricity and magnetism, hydrodynamics and elasticity, these equations 
being deduced finally in terms of a general curvilinear coordinate system. 

In Chapter 7 there is a return to a general Riemannian space, and such topics as relative tensors, 
volume, Green’s theorem and Stokes’ theorem are discussed. Chapter 8 deals with certain non-Rie- 
mannian spaces, in which there is no metric tensor and the tensor derivative is defined through the 


introduction of a linear connection. 


G. E. Hay 

















